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I.   INTRODUCTION 


This  report  is  the  fourth  in  a  series  of  five  volumes  comprising 
the  final  report  for  the  study  entitled  Evaluation  of  Flood  Risk  Factors 
in  the  Design  of  Highway  Stream  Crossings,   authorized  by  the  Federal  High- 
way Administration  (FHWA)  under  Contract  No.  DOT-FH-1 1-7669.  The  overall 
objective  of  the  study  is  to  develop  an  engineering  systems  analysis  tech- 
nique to  enhance  the  decision-making  process  in  the  design  of  highway 
stream  crossings.  The  method  applies  economic  risk  analysis  in  addition 
to  hydraulic  and  hydrologic  factors  to  evaluate  design  alternatives. 

Volume  IV  describes  the  application  of  risk  analysis  to  determine 
the  economically  optimal  design  configuration  for  highway  stream  crossings 
The  analysis  considers  the  sum  of  initial  construction  costs  and  the  flood 
risks,  in  a  range  of  flood  sizes,  for  each  design  alternative.  The  sums 
for  all  design  alternatives  are  compared  and  the  optimal  design  scheme 
selected. 


BACKGROUND 

The  primary  purpose  of  highway  drainage  facilities,  such  as  bridges 
and  culverts,  is  to  move  highway  and  hydraulic  traffic  smoothly,  each  along 
its  own  course  without  interference.  Unlike  highway  traffic,  the  hydraulic 
traffic  cannot  be  controlled  by  enacting  statutes  or  erecting  road  signs. 
Consequently,  conflicts  exist  between  the  highway  system  and  the  waterway. 


The  major  source  of  conflict  arises  from  the  modification  of  natural 
waterways  in  both  horizontal  and  vertical  dimensions  and  the  consequent 
alteration  of  stream  characteristics  at  crossing  sites.  A  major  highway 
system  crossing  a  wide  stream  valley  involves  the  installation  of  bridge 
piers  and  approach  embankments  in  the  stream  waterway.  The  result  of  such 
construction  is  a  change  of  hydraulic  behavior  near  or  at  the  crossing 
site.  Notable  changes  in  hydraulic  behavior  are  (1)  an  increase  in  water 
surface  elevation  upstream  from  the  bridge,  (2)  an  increase  in  flow  velocity 
under  the  bridge,  and  (3)  the  possible  inundation  of  bridge  decks  or 
overtopping  of  approach  embankments.  The  obvious  effects  associated  with 
the  construction  of  highway  culverts  are  ponding  in  the  upstream  areas 
and  increase  in  flow  velocity  downstream  from  the  culvert  exit,  and  in 
some  cases,  roadway  overtopping. 

Depending  upon  the  extent  and  nature  of  stream  cross-section 
modifications  at  the  bridge  site,  a  given  flood  may  or  may  not  result  in 
greater  flood  damage  than  under  preconstruction  conditions.  In  some 
cases  it  is  possible  that  excessive  property  damage  and  heavy  loss  of 
human  life  may  occur  if  the  crossing  is  not  properly  designed.  A  design 
engineer  thus  has  an  extremely  important  task  in  the  process  of  decision- 
making for  his  design.  For  safety  reasons,  he  may  lean  toward  an  over- 
design  alternative,  but  then  find  that  he  has  exceeded  the  available 
project  funds.  Furthermore,  while  it  is  possible  to  predict  the 
probability  of  flood  events  of  various  sizes  over  some  period  of  time, 
it  is  not  possible  to  predict  the  specific  occurrence  of  such  events, 
which  makes  design  decisions  all  the  more  difficult. 

Thus  an  adequate  design  scheme  is  difficult  to  define,  as  it  may 
be  adequate  for  floods  of  magnitude  less  than  the  design  flood,  but  not  for 
those  floods  which  are  larger  than  the  design  flood.  Since  it  is  impracti- 
cal to  design  a  project  for  extreme  flood,  a  certain  amount  of  risk  in  terms 
of  structural  failure,  property  damage  and  loss  of  life  always  exists. 


It  is  the  designer's  objective  to  design  a  project  which  is  least 
costly  to  society  as  a  whole.  To  evaluate  such  costs  the  cost  of  initial 
investment  together  with  the  risk,  or  chance  of  flood  loss,  for  all  flood 
sizes  must  be  considered.  This  report  illustrates  the  procedures  for  per- 
forming such  analysis. 

The  analysis  procedure  developed  in  this  study  applies  only  to 
highway  stream  crossings  on  flood  plains  with  no  constraints.  The  governing 
design  problem  is  assumed  to  be  to  what  extent  the  flow  can  be  constricted 
without  exceeding  some  acceptable  level  of  flood  damage.  The  need  for  such 
analysis  is  evident  from  a  recent  study  by  Chang  (1).  Chang  reports  that 
the  cause  of  bridge  failures  according  to  the  Emergency  Relief  Files  is 
"flow  path  defficiency."  Furthermore,  the  records  of  the  FHWA  indicate 
that  more  than  80  percent  of  the  federal  aid  bridges  built  between  1953 
and  1973  are  over  water.  This  report  describes  a  means  for  assessing  the 
optimum  combination  of  bridge  opening  and  embankment  height  for  a  proposed 
bridge  site.  It  is  recognized  that  there  are  many  bridge  situations  that 
do  not  lend  themselves  to  the  optimization  procedure  developed  in  this 
study,  including  such  cases  as  spans  over  deep  gorges,  or  spans  over  a 
certain  future  meander  cut-off. 


GENERAL  DESIGN  ALTERNATIVES 

Consider  the  case  where  a  highway  crosses  the  flood  plain  of 
a  major  stream.  There  are  several  alternative  schemes  for  passing  the 
traffic  across  the  waterway.  The  most  obvious  is  to  build  a  bridge 
extending  from  bank  to  bank.  The  initial  investment  for  such  a  scheme, 
however,  may  not  be  economically  feasible  if  the  valley  is  two  or  more 
miles  wide. 


Another  alternative  is  to  design  the  highway  stream  crossing  with 
a  combination  of  approach  embankment  and  bridge  if  the  fill  height  is  not 
excessively  great.  While  such  a  combination  results  in  significantly  lower 
initial  construction  costs  (assuming  earthfill  is  less  expensive  per  unit 
length  than  bridge  structure),  the  relationship  of  embankment  height  and 
length  to  bridge  height  and  length  determines  not  only  cost,  but  also  the 
ability  of  the  highway  crossing  to  survive  a  major  flood  undamaged.  Thus 
fill  height  may  vary  considerably  depending  on  the  importance  of  the  highway. 
The  trade-off  is  in  terms  of  lower  initial  cost  versus  some  probability  of 
potential  cost  in  flood  damage  both  at  the  bridge  site  and  upstream.1 

There  are  two  principal  choices  within  this  combination  scheme. 
One  is  to  build  a  high  embankment  to  insure  against  overtopping.  The 
higher  the  embankment,  the  longer  the  bridge  structure  must  be  to  allow 
an  adequate  flood  flow  channel.  This  choice  may  have  a  relatively  high 
initial  cost,  but  a  lower  potential  flood  damage  cost  at  the  site.  The 
second  choice  is  to  build  a  lower  embankment  and  a  shorter  bridge.  While 
this  mode  carries  a  relatively  low  initial  investment,  it  has  a  higher 
potential  cost  in  flood  damage  in  terms  of  roadway  damage  and  traffic 
disruption  due  to  embankment  overtopping.  However,  upstream  backwater 
damage  may  be  lower  under  these  conditions  than  in  the  first  choice  and 
the  damage  to  the  bridge  structure  may  be  reduced  since  flow  velocity  through 
the  bridge  opening  is  lower  in  this  case,  assuming  overtopping. 

Engineering  economics  offers  a  method  to  find  the  most  cost- 
effective  combination  of  embankment  height  and  bridge  length  considering 
the  probability  of  various  flood  magnitudes.  This  method  is  a  decision- 
making tool  for  evaluating  proposed  bridge  construction  investments  in 
dollar  terms. 


Approach  embankments  crossing  flood  plains  hinder  the  flood  flow  passage 
and  have  the  effect  of  raising  upstream  water  levels,  the  so-called 
backwater,  which  may  cause  upstream  property  damage. 


APPLICATION  OF  RISK  ANALYSIS 

On  interstate  projects,  bridges  and  culverts  conventionally  are 
constructed  to  accommodate  floods  of  at  least  a  50-year  frequency  or  the 
greatest  flood  of  record,  whichever  is  greater,  with  runoff  based  on 
land  development  20  years  in  the  future  and  backwater  limited  to  upstream 
property  or  the  highway.  Executive  Order  11296,  issued  in  August  1966, 
and  Flood  Hazard  Evaluation  Guidelines  for  Federal  Executive  Agencies, 
issued  in  May  1972,  require  that  the  elevation  of  the  100-year  flood  at 
project  sites  be  delineated.  Since  the  publication  of  these  guidelines, 
the  FHWA  has  issued  a  new  policy  for  hydraulic  design  of  highway  encroach- 
ments on  flood  plains.  The  new  policy  is  given  in  Appendix  A. 

Both  the  Flood  Hazard  Evaluation  Guidelines  and  the  FHWA  design 
policy  indicate  the  need  for  analyzing  the  risk  of  a  given  project  for 
more  than  a  single  design  flood,  including  the  potential  risks  that  may 
be  derived  from  floods  larger  than  the  design  flood.  Such  considerations 
are  needed  as  is  evident  by  the  recent  severe  effects  of  Storms  Camille 
and  Agnes. 

During  these  two  tropical  storms,  massive  damage  and  inconvenience 
were  reported  in  the  eastern  and  southeastern  United  States.  The  states 
of  Virginia,  Maryland  and  Pennsylvania  suffered  a  total  of  nearly  $75  million 
in  highway  and  bridge  damage  by  Agnes  in  June  1972.  For  this  same  storm 
the  Pennsylvania  Department  of  Transportation  (2)  reported  that  252  bridges 
were  out  of  service  on  that  State's  highway  system,  including  federal -aid 
and  nonfederal -aid  roads.  In  August  1969,  Camille  caused  Virginia  highway 
damage  assessed  at  $19  million,  destroying  or  damaging  133  bridges  and 
closing  25  miles  of  primary  roads.  The  Virginia  Department  of  Highways 
estimates  that  Camille  caused  over  $133  million  to  the  Commonwealth  in 
personal  and  property  damage  (3).  Damages  on  this  scale  warrant  a  closer 
look  at  the  conventional  method  on  which  investment  decisions  for  highway 


stream  crossings  have  been  made.  Logical,  methodical  decision-making 
techniques  must  be  developed  and  implemented  so  that  the  potential  damages 
can  be  weighed  against  initial  cost. 

Risk  factor  analysis  is  a  method  which  has  been  applied  in  other 
fields  for  similar  purposes.  In  risk  analysis,  as  applied  to  stream 
crossing  designs,  the  expected  value  of  flood-related  damages  to  the 
project  itself  and  other  losses  caused  by  the  waterway  construction  are 
evaluated  for  a  complete  spectrum  of  flood  frequencies.  This  value  is 
weighted  and  added  to  the  initial  cost  of  the  project.  With  this  infor- 
mation, an  optimum  design  can  be  determined.  The  basic  theory  of  this 
method  is  to  find  a  design  scheme  which  is  least  costly  to  society 
rather  than  to  design  a  project  for  a  predetermined  flood  occurrence 
interval. 


OBJECTIVE  OF  STUDY 

The  primary  objective  of  this  study  is  to  develop  a  method  based 
on  risk  analysis  to  select  the  most  cost-effective  relationships  among 
the  major  parameters  in  the  design  of  highway  stream  crossings.  This 
method  differs  from  the  conventional  design  method  by  taking  into  account 
flood  damages  to  the  structure,  to  adjacent  properties  from  backwater, 
and  time  loss  of  traffic  delay  during  crossing  overtopping  and  detouring. 
Each  design  scheme  is  evaluated  by  an  annual  total  cost  consisting  of 
annual  cost  of  amortizing  the  structure,  maintenance,  and  the  annual 
expected  value  of  flood  damages  (i.e.,  risk).  The  method  uses  dollar 
costs  as  a  basis  for  comparing  the  results  associated  with  each  set  of 
design  parameters. 

Highway  statistics  indicate  that  approximately  25  percent  of 
U.S.  highway  construction  costs  are  allocated  to  drainage  facilities. 
This  amounts  to  a  cost  ranging  from  $50,000  per  mile  for  secondary  roads 


to  nearly  $1,000,000  per  mile  on  the  interstate  system.  This  high  cost 
warrants  the  development  and  application  of  a  systematic  tool  for 
selecting  the  most  cost-effective  designs  applicable  to  each  site. 


SCOPE  OF  STUDY 

The  basic  concept  of  the  economic/risk  analysis  method  is  to 
integrate  the  hydrologic,  hydraulic  and  economic  factors  to  assess  the 
total  economic  response  for  a  given  design  scheme.  The  annual  total 
economic  response  is  equal  to  the  sum  of  the  annual  initial  cost1   of 
the  project  plus  the  annual  risk.  Risk  is  defined  as  the  expected 
yearly  losses  associated  with  flooding;  this  definition  weights  the 
economic  effects  by  their  probable  outcomes.  Flood  losses  include  dam- 
age incurred  to  bridge  and  embankment  structures,  backwater  damage  to 
property  adjacent  to  the  highway  crossing,  loss  of  time  to  highway  users 
when  the  bridge  is  impassable  during  embankment  overtopping  and  the  sub- 
sequent repair  period,  and  additional  cost  due  to  detouring. 

All  the  analyses  are  to  be  carried  out  on  a  digital  computer 
in  FORTRAN.  While  the  method  of  analysis  is  intended  to  be  as  compre- 
hensive as  possible,  certain  assumptions  have  been  made  in  the  model 
development.  In  general,  the  method  is  defined  within  the  following 
framework: 

1.  The  computation  of  water  surface  elevation  and  flow 
velocities  throughout  the  flow  region  is  performed 
by  a  two-dimensional  Finite  Element  Model.  These 
hydraulic  computations  give  solutions  for  steady 
state  peak  flow  conditions. 

2.  Flood  hydrographs  and  their  return  periods  are  the 
hydrologic  data  needed  to  conduct  the  analysis. 

2Defined  on  page  14  as  the  total  construction  cost  spread  over  the 
amortization  period  of  the  structure. 


3.  Bridge  damage  and/or  loss  by  floating  debris  and  flow 
inundation  are  considered  in  the  model. 

4.  Property  damage  due  to  backwater  resulting  from  the 
bridge  construction  is  represented  by  a  stage-damage 
curve,  as  is  conventional. 

5.  The  method  of  analysis  presented  in  this  report  cannot 
be  applied  to  assess  the  case  of  irreducibles,  such  as 
defense  or  emergency  evacuation  routes  or  for  the  case 
where  severe  budget  constraints  are  imposed  on  design. 

6.  Economic  factors  must  be  estimated  on  a  case-by-case 
basis. 


II.  SITE  DATA  COLLECTION 


The  first  step  of  an  engineering  analysis  for  a  proposed  high- 
way stream  crossing  is  to  collect  the  necessary  data  for  that 
bridge  site  .  Data  to  be  gathered  for  each  site  are: 


1.  Location  map  showing  the  proposed  highway  alignment, 
embankment  and  structures,  reach  of  river  to  be  affected 
by  embankment  encroachment,  and  existing  embankment 
encroachment  and  highway  structures,  if  any. 


Vicinity  map  showing  flood  flow  patterns,  cross 
sections  of  stream,  location  of  proposed  bridge 
opening,  alignment  of  piers,  skew  of  crossing,  bends 
and  stream  meanders,  and  vegetation  types  and  density 
on  flood  plains. 


Description  of  existing  bridge  structures  both  upstream 
and  downstream  from  crossing  or  encroachment.  This 
should  include: 

a.  Type  of  bridge,  including  span  lengths  and  pier 
orientation. 

b.  Foundation  type  such  as  spread  footing  or  piling; 
foundation  depth. 

c.  Scour  history  at  abutments;  stream  aggradation, 
degradation. 

d.  Stream  cross  section  beneath  structures,  noting 
stream  clearance  to  superstructure  and  skew  with 
direction  of  current  during  floods. 

e.  All  available  flood  history  information,  including 
highwater  marks  with  dates  and  elevation,  nature 
of  flooding,  damages  and  source  of  information. 

f.  Photographs  showing  existing  structures,  past 
floods,  main  channels  and  flood  plains. 


g.  Information  on  nature  of  drift,  ice,  streambed,  bank 
stability,  and  land  use. 


4.  List  of  factors  affecting  water  stage  at  bridge  site. 

a.  Highwater  from  other  streams. 

b.  Reservoirs  -  existing  or  proposed  and  approximate 
date  of  construction. 

c.  Flood  control  projects  (give  status). 

d.  Tides. 

e.  Other  controls. 

5.  Geologic  data  at  proposed  bridge  site. 

6.  Hydrologic  data. 

a.  Drainage  area  above  the  proposed  bridge  crossing  site. 

b.  Streamflow  records  at  the  site  and/or  at  nearby 
stations. 

c.  Characteristics  of  drainage  basin  including  basin 
slope,  soil  type,  and  rainfall. 

7.  Traffic  data. 

a.  Average  daily  traffic  (ADT) . 

b.  Vehicular  diversity  such  as  fraction  of  ADT  in 
passenger  cars,  commercial  vehicles,  etc. 

c.  Number  of  occupants  per  vehicle. 

d.  Speed. 

e.  Available  detour  routes  and  length  of  detour. 

f.  Accident  data. 

8.  Economic  data. 

a.  Cost  of  labor  and  construction  materials. 

b.  Average  personal  income. 

c.  Type  of  land  development  alongside  the  stream,  including 
agricultural,  manufacturing,  residential,  retail 
business,  selective  services  and  wholesale  business. 
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There  are  many  possible  sources  for  the  data  listed  in  the  preceding 
paragraphs.  Depending  upon  the  type  needed,  data  may  be  collected  by 
federal,  state  or  private  agencies.  Typically,  data  types  1  and  2  are 
obtainable  from  the  U.  S.  Geological  Survey  (USGS)  guadrangle  maps. 
To  allow  a  more  detailed  description  of  the  site,  aerial  photographs  and 
specially  prepared  maps  with  1-  or  2-foot  contours  must  be  used.  In  some 
cases,  cross  sections  normal  to  flood  flow  are  used  in  lieu  of  topographic 
maps.  As  a  minimum,  three  sections  are  required:  one  upstream,  one  at  the 
site  crossing,  and  one  downstream.  The  type  and  distribution  of  ground 
cover  on  the  flood  plains  may  be  obtained  from  aerial  photos  and  site 
visits. 

Geological  data  are  generally  obtainable  from  either  USGS  and/or 
state  geological  surveys  or  departments.  Boring  may  be  required  if  no 
geological  data  can  be  obtained  from  these  or  other  sources. 

Traffic  data  are,  of  course,  gathered  by  the  state  highway  depart- 
ments. Most  of  the  economic  data  are  obtainable  from  national  census  data. 

Streamflow  data  have  been  collected  at  gaging  stations  throughout 
the  country  and  published  by  the  USGS.  The  U.  S.  Weather  Service  provides 
rainfall  data  for  various  locations  throughout  the  country. 

At  a  proposed  bridge  site  the  streamflow  record  frequently  is 
either  nonexistent  or  for  only  a  limited  period.  With  such  limited 
information  it  is  difficult  to  obtain  accurate  hydro! ogic  data  for 
use  in  analysis.  In  order  to  overcome  this  shortcoming  the  FHWA  is 
currently  sponsoring  research  projects  to  facilitate  better  prediction 
computations  methods  for  peak  flow  from  small  watersheds.  It  is  anticipated 
that  once  these  results  become  available  they  will  be  incorporated  into 
the  model  developed  in  this  report.  For  the  purpose  of  model  development 
in  this  study,  however,  the  hydrologic  data  which  are  available  at  the 
present  time  have  been  used  for  analysis. 
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III.  METHOD  OF  ANALYSIS 


The  general  concept  and  approach  for  applying  risk  analysis  to 
a  series  of  design  schemes  for  a  highway  stream  crossing  includes  the 
computation  of  annual  construction  cost  and  expected  yearly  risk  asso- 
ciated with  floods  for  selected  design  schemes.  The  sum  of  the  annual 
construction  cost  and  the  expected  flood-related  loss,  or  risk,  is  the 
economic  response  or  expected  tangible  total  costs  of  the  stream  crossing 
to  the  highway  users. 

Figure  1  shows  the  overall  logic  sequence  and  analysis  procedure 
for  evaluating  each  design  scheme.  For  each  scheme  the  procedure  involves 
five  major  steps: 

1.  Calculate  annual  construction  costs, 

2.  Perform  hydraulic  computation, 

3.  Estimate  embankment  erosion  and  scour  under  bridge, 

4.  Compute  losses  associated  with  structural  and  property 
damages,  and  losses  incurred  from  traffic  delay,  and 

5.  Weight  losses  with  flood  probabilities  to  determine  risk. 

The  construction  costs  are  computed  for  each  design  scheme  by  using 
either  a  quantity-times-unit-cost  method  or  empirical  methods.  Structural 
damage  is  assumed  to  be  directly  related  to  the  extent  of  embankment  erosion 
caused  by  flow  overtopping,  degree  of  debris  clogging  at  the  bridge  opening, 
depth  of  bridge  inundation,  scour  caused  by  the  excessive  velocity  under  the 
bridge  opening,  and  the  washout  of  roadway  from  overtopping.  Flood  back- 
water damage  to  property  is  derived  from  a  flow-  or  stage-damage  function 
formulated  for  each  crossing  site.  Traffic-related  losses  include  the 
cost  of  lost  time,  increased  running  costs,  accident  losses  on  the  detour, 
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Figure  1.  Logic  Flow  of  Risk  Analysis 
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and  accident  losses  due  to  an  unexpected  obstacle  or  barricade  placed  at 
the  stream  crossing  site  when  a  failure  of  the  roadway  occurs. 

The  estimated  structural  damage,  flood  damage,  and  traffic-related 
losses  are  multiplied  by  the  probability  of  yearly  flood  occurrence  to 
obtain  the  risk  for  each  flood  event.  Summing  the  risks  over  the  set  of 
flood  events  selected  for  study  yields  the  risk  component  of  the  total 
economic  response. 


ANNUAL  CONSTRUCTION  COST 

The  construction  costs  for  a  highway  crossing  normally  consist  of: 

Bridge  superstructure, 

Bridge  substructure  and  excavation, 

Approach  embankments, 

Roadway  pavement,  and 

Protective  measures  (spur  dikes,  riprap). 

In  general,  the  total  construction  costs  are  governed  by  the  type  of 
bridge,  bridge  length,  bridge  width,  clearance,  span,  geology  of  the 
river  bed,  height  of  the  embankment,  type  of  highway  (e.g.,  primary  or 
secondary),  and  width  of  the  valley  at  the  crossing  site.  In  the  model 
construction  costs  are  calculated  as  the  sum  of  the  following  cost  components 

(1)  Roadway  pavement, 

(2)  Embankment, 

(3)  Bridge,  including  the  superstructure,  substructure  and 
excavation,  and 

(4)  Protective  measures. 
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1.  COST  OF  ROADWAY  PAVEMENT 

An  estimate  of  the  cost  of  this  parameter  is  calculated  by  multi- 
plying the  total  length  of  the  roadway  between  the  beginning  and  end 
stations  of  the^valley  crossing,  excluding  the  bridge  length,  by  the  cost 
per  linear  foot  ($/L-ft)  of  roadway.  It  is  assumed  that  (1)  the  cost 
factor  includes  materials  and  labor,  and  (2)  the  cost  of  roadway  pavement 
over  the  bridge  is  included  in  the  bridge  cost. 

2.  COST  OF  EMBANKMENT 

Calculation  of  construction  cost  for  approach  embankments  is  made 
on  a  unit  volume  basis.  Given  grade  and  roadway  elevation  along  the 
highway  center  line,  the  embankment  top  width,  and  upstream  and  downstream 
side  slopes,  the  quantity  of  earth  fill  can  be  readily  computed.  The 
total  cost  is  computed  by  multiplying  the  quantity  of  fill  volume  by  the 
unit  cost  ($/cu  ft).  Figure  2  shows  the  geometry  of  a  typical  river 
crossing. 

In  the  model  the  volume  of  the  i   element  between  stations  i  and 
i+1  is  calculated  by 

Vol.  =  rw  VLHi±L  L.  +l(sl+S2)(H?  +  H|+1)Li    (1) 

where   H. ,  H.+1  =  fill  height  at  stations  i  and  i+1,  respectively, 

RW  =  width  of  roadway, 
L.  =  element  length,  and 


S-, ,  Sp  =  side  slope  of  fill 
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3.     BRIDGE  COSTS 

The  calculation  of  the  total  Risk  associated  with  a  structure  requires 
an  accurate  evaluation  of  the  cost  of  construction  for  the  structure.  Bridge 
costs  can  be  divided  into  the  cost  of  the  following  major  components:  bridge 
deck,  piers,  abutments,  footings  and  piles.  Each  of  the  components  can  be 
expressed  either  as  a  function  of  the  bridge  length,  embankment  height,  expected 
scour  depth,  as  in  the  case  of  piers,  or  a  constant  cost.  For  example,  bridge 
deck  cost  is  related  to  the  length  of  the  span  while  pier  cost  can  be  expressed 
as  a  function  of  expected  scour  depth.  Conceptually,  a  bridge  costing  function 
can  take  the  form: 

Cb  "  Csup  +  Csub  ■  «L"H'  <*) 

where 

C.  =  total  costs 

Csup  =  superstructure  costs 

C  .  =  substructure  costs 

L  =  bridge  length 
H  =  fill  height 

A  portion  of  the  substructure  costs  includes  the  costs  due  to  the  depth 
piers  must  be  driven  in  anticipation  of  expected  scour.  Presently,  the  model 
allows  one  pier  depth  to  be  input  for  all  the  model  runs  using  the  scour  option. 
If  a  conservative  approach  to  the  analysis  is  followed  this  requires  this  depth 
and  its  associated  costs  to  be  calculated  for  the  narrowest  opening  and  maximum 
flow  and  embankment  elevation. 
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This  is  suggested  for  two  reasons: 

1.  Bridges  are  not  purposely  designed  to  fail  under  any 
expected  conditions  and  to  allow  the  structure  to  fail 
during  analysis  represents  a  condition  not  followed  in 
standard  practice. 

2.  Since  only  one  pier  depth  is  permitted,  a  conservative 
approach  to  the  structural  design  of  the  bridge  requires 
the  piers  be  designed  for  the  worst  case,  which  is  maximum 
flow,  maximum  embankment  elevation,  and  minimum  opening 
size. 

Ideally  a  third  factor  besides  embankment  height  and  bridge  opening  in 
the  analysis  would  be  pier  length  and  its  costs.  However,  this  considerably 
increases  the  complexity  of  the  problem.  Useful  information  could  be  obtained 
by  a  more  detailed  study  of  this  problem  and  can  easily  be  incorporated  into  the 
Risk  analysis  at  some  future  time. 

The  costing  method  presently  employed  allows  the  user  two  options.  He 
can  either  input  a  cost  function  for  each  bridge  design  or  allow  the  model  to 
use  its  own  internal  functions  as  described  in  Table  1.  Option  1  is  formulated 
below: 

Cb  =  f(L)  +  G*d1  +  K  (3) 

where 

f(L)  =  B]  +  B2  *  L 

and  K  =  B3  +  B. 

B.|  =  constant  coefficient  for  bridge  superstructure  cost 
B2  =  coefficient  of  bridge  length 
B3  =  pier,  footing  and  abutment  costs 
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B*  =  spur  dike  costs 

G  =  cost/ft  of  pile  (optional) 
d,  =  length  of  pile  (optional) 

Equation  3  constitutes  the  basic  bridge  cost  equation  for  the  model. 
Due  to  the  variation  in  bridge  type  and  the  variation  in  costs  of  construction 
materials  and  labor  from  location  to  location,  it  is  generally  difficult  to 
formulate  an  equation  of  C.  sufficiently  convenient  for  modeling  use.  Conse- 
quently, the  user  must  provide  the  coefficients  as  inputs  to  the  model  for 
determination  of  bridge  cost.  It  should  be  noted  that  Equation  3  also  provides 
the  user  with  the  opportunity  to  input  the  total  bridge  cost  directly,  if  such 
information  is  available.  This  may  be  accomplished  by  setting  the  coefficients 
B«  through  B«  to  zero  and  inputing  B-,  as  the  known  total  bridge  cost. 

« 

A  second  option  is  available  for  those  users  who  are  only  interested 
in  approximate  bridge  cost  estimation  or  who  are  unfamiliar  with  bridge  cost 
computations.  The  second  option,  described  below  is  permanently  included  in 
the  model-.  The  information  provided  in  the  second  option  was  derived  based  on 
statistical  data  of  bridge  costs  throughout  the  United  States  collected  over  a 
period  of  21  years.  The  bridge  cost  equation  used  in  the  routine  was  developed 
during  the  course  of  this  study  based  on  a  regression  analysis  of  costs  on 
bridge  lengths.  Data  were  provided  by  the  FHWA  (4),  which  has  collected 
extensive  files  of  cost  data  through  the  federal  aid  highway  programs.  The 
bridge  cost  data  contained  in  the  FHWA  file  (4)  for  the  period  1953-1972  were 
regressed  in  the  form  of 


where 


Cb  =  aQ  +  alXl  (4) 


a  =  constant 
o 

a,  =  coefficient  for  length  (x-,) 
X-.  =  bridge  length 
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The  results  are  shown  in  Table  1. 

This  method  of  estimating  cost  is  not  accurate  for  bridges  of  less  tfa&n 
0.1  mile  in  length  since  scatter  observed  in  the  cost  data  was  too  large. 

Table  1.  Bridge  Cost  Regressed  on  Bridge  Length1 

Length  an       a.      a 

(Miles)  o        4       yx         Y 


.004  -  .009 

750. 

-1421. 

2611. 

.12 

.01  -  .019 

1165. 

-295. 

1295. 

.05 

.02  -  .039 

9752. 

1593. 

687. 

.46 

.04  -  .099 

1467. 

1722. 

665. 

.51 

.10  -  .199 

4830. 

1547. 

163. 

.91 

>.20 

5080. 

1608. 

184. 

.89 

where  a   =  standard  error,  and 

y  =  corre Tation  coefficient. 

1Costs  are  based  on  1972  price  level.  Adjustment  for  inflation  has  been  con- 
sidered for  the  period  studied. 

4.     PROTECTIVE  MEASURES 


Constriction  of  flood  plain  by  approach  embankment  invariably 
induces  a  lateral  movement  of  flood  water  from  the  embankment  face 
toward  the  bridge  opening.  As  this  flow  joins  the  main  flow  a  high  flow 
concentration  in  the  vicinity  of  the  abutment  occurs.  This  high  flow 
concentration  may  produce   violent  turbulence  as  it  enters  the  constric- 
tion, thus  creating  scour  at  the  abutment  toe. 

An  effective  means  of  preventing  abutment  scour  is  the  installa- 
tion of  spur  dikes  attached  to  the  abutments.  The  length  of  spur  dikes 
is  generally  governed  by  the  fraction  of  flow  on  the  flood  plain  and  the 
velocity  of  the  main  flow.  A  method  for  determining  the  length  of  spur 
dikes  is  presented  in  Reference  5.  The  height  of  spur  dikes  is  determined 
based  on  water  depth  in  the  flood  plain. 


20 


The  side  slope  of  spur  dikes  is  governed  by  the  type  of  material 
to  be  used  for  spur  dike  construction.  Once  the  geometry  of  the  dikes 
is  known,  the  quantity  of  material  and  hence  the  cost  can  be  obtained. 
This  cost  is  added  directly  to  the  construction  cost  of  xtne  bridge 
under  option  1  of  calculating  bridge  cost  and  is  included  in  the 
regression  equation  used  in  option  2. 

At  some  bridge  sites  the  use  of  riprap  is  required  to  protect 
the  embankments  and  abutments  from  erosion.  The  method  of  computing 
the  size  of  riprap  materials  is  given  in  Reference  6.  The  cost  of 
the  riprap  is  added  to  the  construction  cost  of  the  bridge  directly 
under  option  1  of  calculating  bridge  cost  and  is  included  in  the  regres- 
sion analysis  employed  under  option  2. 

The  four  major  construction  cost  components  (roadway  pavement, 
embankment,  bridge  superstructure,  substructure  and  excavation,  and 
protective  measures)  are  summed  to  determine  the  total  initial  construc- 
tion cost  of  the  bridge  design.  The  annual  construction  cost  is  obtained 
by  multiplying  the  total  initial  construction  cost  by  the  capital  recovery 
factor  for  the  amortization  period  (life  of  the  structure)  and  the  appro- 
priate interest  rate.  The  capital  recovery  factor  is  expressed  as: 

CRF \ +  i  (5) 

(1  +  i)n  -  1 


i 


where 


i  =  appropriate  interest  rate  and 
n  =  amortization  period  in  years. 
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HYDROLOGY 

A  detailed  description  of  the  hydro! ogic  aspects  of  highway  stream 
crossing  sites  is  beyond  the  scope  of  this  study.  It  is  assumed  that 
hydro! ogic  data  necessary  for  risk  analysis  are  available  for  study  sites 
of  interest. 

Two  types  of  hydrologic  data  are  essentia!  for  risk  analysis.  The 
first  type  is  the  flood-frequency  curve,  and  the  second  type  is  flow 
hydrographs  of  each  individual  flood  event.  The  flood-frequency  curve 
provides  data  on  peak  discharges  and  the  associated  probability  of 
occurrence  needed  for  computing  flood  risks. 

A  waterway  constriction  controls  the  flow  with  an  effect  similar 
to  that  of  a  reservoir  spillway  or  sluice  gate.  The  stage-discharge 
relationship  is  altered  once  the  control  structure  is  in  place.  Dynamic 
routing  of  flows  before  and  after  this  constriction  shows  the  difference 
in  hydrographic  response.  However,  as  dynamic  routing  is  not  within  the 
scope  of  this  study,  inflow  hydrographs  and  the  stage-discharge  curve 
at  the  bridge  site,  with  bridge  in  place,  are  used  to  estimate  the 
duration  of  embankment  overtopping. 

The  runoff  hydrographs  for  this  study  are  assumed  to  be  triangular 
in  shape.  In  order  to  define  these  hydrographs,  three  parameters  are 
required:   (1)  time  to  peak,  T  ;  (2)  flood  duration,  T.  ;  and  (3)  the 
peak  flow,  Q  .  A  general  representation  of  the  assumed  runoff  hydrograph 
is  shown  in  Figure  3. 

Time  to  peak  is  estimated  using  the  drainage  basin  characteristics 
incorporated  in  a  design  method  used  by  the  U.  S.  Bureau  of  Reclamation  (7) 
Time  to  peak  is  approximated  by  the  time  of  concentration,  Tc,  defined 
as  the  travel  time  of  the  runoff  from  the  hydraulically  most  distant 
point  to  the  bridge  site.  The  empirical  method  presented  by  the  Soil 
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Figure  3.  General  Shape  of  Assumed  Hydrograph 
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Conservation  Service  (8)  for  determining  the  time  of  concentration  is 


where 

L  =  length  of  longest  watercourse  in  miles, 

H  =  elevation  difference  in  feet,  and 

T  =  time  of  concentration  in  hours. 

The  Bureau  of  Reclamation  method  is  used  to  determine  the  length  of  time 

between  the  peak  flow  and  the  end  of  the  hydrograph.  This  time  period 

is  designated  T  in  Figure  3.  For  a  given  watershed,  the  relationship 
between  the  time  to  peak  (T  )  and  Tr  is 

r 


f=K 


The  value  of  the  constant,  K,  is  computed  from  recorded  hydrographs 
Analyses  by  the  Soil  Conservation  Service  (8)  suggest  1.67  as  an  average 
value  for  this  constant  for  use  on  ungaged  watersheds.  This  constant 
substituted  into  the  above  equation  produces  Equation  7,  which  is  used 
in  the  RISK  module  to  determine  T  . 


Tr  =  1.67  Tp  (7) 


In  summary,  we  have  the  following  relationships: 

t   t   /n.9  L3\  -385 

'd  "  'C  "      u 


Tr  =  Tb  -  Tp  =  1.67  Tp 
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The  final  parameter  required  to  determine  the  hydrograph  is  the 
peak  flow  (QD).  The  range  of  peak  flows  associated  with  different  recur- 
rence  intervals  is  determined  by  fitting  runoff  data  at  the  site  to  a 
Gumbel  plot.  A  Gumbel  plot  is  a  linearized  graph  of  relative  flood  peak 
magnitudes  versus  the  average  recurrence  interval  in  years.  The  axes  are 
divided  in  such  a  manner  that  the  storm  frequency  distributions  plot  as 
an  approximately  straight  line.  Figure  4  shows  a  Gumbel  plot  for  the 
Tallahalla,  Mississippi,  bridge  site.  The  data  for  the  figure  were 
recorded  at  Laurel,  Mississippi,  the  nearest  station  to  the  Tallahalla 
bridge  site.  For  demonstration  purposes,  the  data  were  not  corrected 
for  drainage  area  difference  since  the  station  at  Laurel  is  within  five 
miles  of  the  bridge  site,  with  no  major  tributaries  draining  into  the 
five-mile  reach. 


HYDRAULIC  COMPUTATIONS 

As  mentioned  previously  the  construction  of  bridge  approach 
embankments  onto  flood  plains  results  in  the  constriction  of  flood  water- 
ways. This  could  cause: 

1.  Increase  in  water  surface  elevation  over  normal  stage 
Ci.e.,  backwater)  on  flood  plains  and  adjacent  lands 
near  the  bridge  crossing  site. 

2.  Overtopping  of  approach  embankments  during  flooding. 

3.  Increase  in  flow  velocity  through  the  bridge  opening. 

All  of  the  above  three  cases  could  produce  varying  degrees  of 
flood  loss.  Backwater  flood  damages  are  incurred  to  real  properties 
such  as  commercial,  industrial  and  residential  developments.  Embank- 
ment overtopping  causes  losses  ranging  from  traffic  delay  to  structural 
failure.  An  excessive  increase  in  velocity  through  the  bridge  opening  up- 
sets the  hydraulic  regime,  eroding  the  stream  bed  and  causing  local  scour 
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around  bridge  piers.  In  order  to  assess  the  damage  that  a  given  flood 
could  produce,  the  hydraulic  information  associated  with  that  flood  must 
be  properly  computed  under  a  preselected  set  of  bridge  design  schemes. 

In  the  real  world  flood  flow  is  a  dynamic  phenomenon.  The 
variation  in  flood  stage  and  velocity  components  is  both  temporal  and 
spatial.  Unfortunately,  techniques  for  dynamically  routing  multi- 
dimensional flow  are  still  quite  limited  in  both  applicability  and  practi- 
cability. In  the  case  of  bridge  hydraulics  the  ultimate  interest  of  the 
designer  lies  in  the  peak  flow  conditions.  For  this  reason  the  hydraulic 
computation  is  necessary  for  steady  state  peak  flow  conditions  only. 

Data  required  for  the  hydraulic  computation  are: 

1.  Topography  or  cross  section  of  stream  near  bridge  site, 

2.  Aerial  photo, 

3.  Vegetation  map, 

4.  Stage-discharge  diagram, 

5.  Flood-frequency  curve,  and 

6.  Bridge  configuration—opening  size,  skewness,  eccentricity. 

The  hydraulic  description  is  obtained  by  solving  the  set  of  govern- 
ing partial  differential  equations  for  the  flow  region,  assuming  prescribed 
boundary  conditions.  The  set  of  partial  differential  equations  consists 
of  the  equations  of  motion  and  continuity  in  a  two-dimensional  flow  region. 
A  numerical  technique,  the  Finite  Element  Method,  is  applied  to  solve 
the  governing  flow  equations.  Throughout  this  report  the  hydraulic  model 
is  referred  to  as  the  Finite  Element  Model  (FEM).  The  hydraulic  output 
of  FEM  is: 


1.  Water  surface  elevation  at  various  locations  throughout 
the  flow  regime, 

2.  Velocity  components  at  various  locations  along  the 
direction  of  the  stream  and  the  transverse  direction, 
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3.  Velocity  through  the  bridge  opening,  and 

4.  Flow  overtopping,  if  any,  of  the  approach  embankments. 

The  hydraulic  computation  is  carried  out  for  each  flood  of 
interest.  The  hydraulic  data  for  each  flood  are  stored  on  hydraulic 
tapes  for  the  computation  of  flood  risk. 

The  Finite  Element  Model  is  described  in  detail  in  Volume  III 
of  this  series.  Chapter  VI  of  that  volume  describes  the  results  of 
the  hydraulic  solution  for  the  Tall  aha! la  Creek  bridge  site. 


CALCULATION  OF  LOSSES 


STRUCTURAL  DAMAGES 


At  a  highway  stream  crossing  site,  flood-related  structural 
damage  results  from  one  or  more  of  the  following: 

1.  Damage  to  bridge  superstructure  due  to  debris  and 
inundation  of  bridge  deck. 

2.  Embankment  erosion  due  to  overtopping. 

3.  Scour  under  bridge  piers  and  abutments. 

1.     Damage  to  Bridge  Superstructure 

When  the  bridge  deck  is  inundated  by  flood  water,  the  super- 
structure is  subject  to  two  types  of  force,  in  addition  to  its  own 
weight.  In  the  horizontal  direction  the  superstructure  is  subject 
to  a  thrust  caused  by  the  dynamic  forces  of  flow  acting  on  the  upstream 
face  of  the  bridge.  In  the  vertical  direction  the  bridge  is  acted  upon 
by  a  buoyancy  force.  Because  of  the  buoyancy  the  effective  weight  of 
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a  concrete  bridge  is  reduced  to  about  60  percent  of  its  weight  in  air. 
The  combination  of  fluid  dynamic  forces  and  buoyancy  may  push  or  lift 
the  superstructure  from  the  abutments  and  piers,  causing  bridge  failure. 

There  are  other  factors  which  may  aggravate  the  bridge  failure: 

•  First,  the  accumulation  of  trash  at  the  bridge  causing 
changes  in  the  flow  patterns  and  thus  unbalancing' the 
forces  acting  on  the  bridge; 

•  Second,  the  impact  of  large  floating  debris  striking 
the  bridge;  and 

t  Third,  the  entrapment  of  air  under  the  deck  between  girders 
resulting  in  further  reduction  of  the  effective  weight  of 
the  bridge  superstructure. 

It  is  generally  difficult  to  predict  the  effects  of  flood 
forces  on  bridge  superstructure  and  the  extent  of  bridge  damage  caused 
by  superstructure  inundation,  particularly  since  data  on  the  yield  and 
nature  of  flood  debris  are  not  easy  to  obtain.  Although  bridge  failures 
due  to  flooding  of  bridge  decks  have  been  reported  during  major  floods, 
this  information  is  insufficient  for  modeling  purposes.  Further 
research  is  urgently  needed  in  this  respect. 

While  bridge  damage  from  inundation  of  superstructure  is  diffi- 
cult to  assess,  it  is  too  important  to  be  ignored  in  assessing  flood 
losses.  For  lack  of  information  the  present  model  uses  the  following 
equation  to  compute  the  bridge  damage,  L  : 

k  =  a  (8) 


where 


C  =  a  coefficient  and 

y  =  submergence  of  bridge  deck. 
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Values  of  t,   may  be  obtained  for  specific  bridge  sites  from  past 
flood  records.  If  c  is  not  specified  in  the  model  input,  it  will  be 
set  to  zero  by  the  program.  This  portion  of  the  algorithm  may  be  modi- 
fied when  a  more  refined  method  for  assessing  bridge  damage  becomes 
available. 

2.     Embankment  Erosion 

Embankment  erosion  caused  by  flow  overtopping  has  not  been 
studied  on  a  broad  scale.  This  is  surprising  considering  the  number  of 
earthfill  dams  and  dikes  which  have  the  potential  for  producing  catastro- 
phes if  they  fail.  While  the  washout  of  a  bridge  approach  embankment 
will  probably  not  result  in  catastrophic  damage,  it  will  cause  traffic 
interruption  and  bridge  structural  damage.  In  order  to  assess  these 
losses,  the  mechanism  of  embankment  erosion  during  overtopping  must  be 
known.  Unfortunately,  existing  literature  can  provide  only  limited 
information  on  the  subject.  A  brief  description  of  related  work  on 
embankment  erosion  follows. 


•  A  series  of  model  tests  by  Posey  (9)  to  study  erosion 
on  uniform  and  graded  riprap  layers  for  various  flow 
conditions;  the  tests  were  inconclusive. 

•  Laboratory  and  field  experiments  by  Tinney  and  Hsu  (10) 
to  demonstrate  the  feasibility  of  using  a  fuse  plug  in 
the  spillway  of  a  major  dam.  In  this  project,  a  pilot 
channel,  ten  feet  wide  and  11.5  feet  high  with  an  invert 
three  feet  lower  than  the  crest,  starts  the  washout  of 
the  fuse  plug.  In  a  1:2  scale  model,  the  breach  time  for 
the  pilot  channel  was  one  minute  27  seconds.  The  lateral 
erosion  rate  was  5.6  feet  per  minute  for  uniform  material 
and  1.4  feet  per  minute  for  well  graded  material.  A  uni- 
form, cohesive  rock  material  was  used  at  the  pilot  channel. 

t  A  method  developed  by  Cristofano  (11)  for  computing  erosion 
rate  in  an  earth  dam  embankment  failure,  assuming  that 
overflow  through  a  shallow  notch  of  a  given  width  had 
started.  The  rate  of  additional  deepening  of  the  notch  or 
erosion  of  fill  material  is  computed  by 
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SJ551L-  «  K  e"x  (9) 

^water 

where 

Q  .-j  =  volume  of  soil  eroded  in  each  time 
period, 

Q. ,a+-^  =  volume  of  water  discharged  each 
water   time  period, 

K  =  constant  of  proportionality, 
1  in  this  case, 

e  =  base  of  natural  logarithmic  system,  and 

x  =  b/H  tan  e 

where   b  =  base  length  of  overflow 

channel  at  any  given  time, 

H  =  hydraulic  head  at  any  given 
time,  and 

0  =  developed  angle  of  friction 
of  soil . 

•  Application  of  Equation  9  by  Newton  and  Cripe  (12)  to  analyze 
the  embankment  breaching  in  safety  studies  of  Tennessee  Valley 
Authority  (TVA)  nuclear  power  plants. 

In  the  absence  of  a  satisfactory  erosion  equation,  the  following 
empirical  equation  developed  in  an  earlier  WRE  study  (13)  is  used  in 
this  study: 

E  =  aV6  (.10) 

where 

E  =  erosion  in  tons/ft/day, 
V  =  velocity  in  ft/sec,  and 
a,  6  =  empirical  constants. 

Estimates  of  a  and  3  are  .25  and  3.8j  respectively.  These  values  repre- 
sent a  compromise  of  the  two  extremes  of  a  cohesive  and  cohesionless  soil, 
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Figure  5  illustrates  the  erosion  mechanism  assumed,  i.e.,  the  eddies  of 
flow  turbulence  and  bed  shear  on  the  embankment  surface.  Once  an  embank- 
ment is  overtopped,  the  head  differential  between  the  upstream  and  down- 
stream sides  of  the  embankment  is  generally  small. 

The  velocity  in  Equation  10  is  based  on  the  assumption  that 
critical  flow  exists  over  the  roadway.  In  view  of  the  small  difference 
in  head  across  the  embankment,  this  assumption  may  be  conservative.  The 
roadway  is  assumed  to  act  as  a  broad  crested  weir  in  the  calculation  of 
the  critical  velocity.  For  bridges,  V  is  assumed  to  be  the  velocity  of 
the  water  as  it  crosses  a  broad  crested  weir;  this  is  the  critical 
velocity  which  corresponds  to  one-third  of  the  potential  energy  associated 
with  the  depth  of  the  water  at  the  weir.  This  bridge  assumption  recognizes 
the  small  water  surface  differences  above  and  below  the  structure  under 
conditions  of  complete  flooding. 

This  embankment  erosion  formulation  assumes  that  after  an  erosion 
threshhold  velocity  (V  )  is  reached,  embankment  erosion  can  be  represented 
by  a  triangular  approximation,  with  area  A  .  Values  for  V  for  various 
soils  are  given  in  Figure  6.  Given  the  area  eroded,  the  amount  of  damage 
to  the  roadway  can  be  calculated  and  a  percent  of  the  total  embankment 
erosion  can  be  calculated  as  follows:  Given  a  section  of  embankment 
length  L,  the  total  volume  of  the  section  is 

V.  =  (B  +  S-H)  H,  L  (11) 


where 


B  =  embankment  width,  and 
H  =  fill  height. 
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\< —  B  — >\ 


v  =  /l 


gd 


=  A  =  [    a 

w  J 


Vp  L  dt 


0 


where  T  =  time  of  overtopping, 
L  =  roadway  length,  and 
d  =  head  above  roadway. 


Figure  5.  Erosion  Mechanism 
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From  Equation  10  the  total  volume  of  erosion  during  the  overtopping  event 
is  calculated  as 

E  =  /  aV3Ldt 
T 

=  aV3LT  =  |  H  L  (12) 


where       T  =  the  time  of  overtopping,  and 
b  =  eroded  roadway  width. 

Thus  the  percent  of  embankment  erosion  is 

Pj  -7  *  100  (13) 

and  the  percent  of  roadway  washout  is 

P2  „bxl0O.21  n4) 

Reiterating,  damage  to  the  approach  embankments  is  of  two  types: 

1.  Erosion  of  the  fill,  and 

2.  Damage  to  the  roadway. 

Since  bridge  damage  is  considered  independently  of  roadway  and  fill  damage, 
this  analysis  does  not  include  abutment  undermining  by  erosion  and  the 
resultant  loss  of  particular  spans.  It  is  assumed  that  bridges  constructed 
in  a  flood  plain  are  designed  to  use  approach  embankments  in  much  the  same 
way  as  an  overflow  spillway  is  used  in  dam  construction. 
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The  economic  loss  incurred  by  damage  to  the  site  (Ls)  is  computed  as: 

Ls  =  (PjCf  +  P2Cr)  Ca  (15) 

where       Cf  =  original  cost  of  fill, 

Cr  =  original  cost  of  roadway,  and 
Ca  =  cost  adjustment  factor. 

The  cost  adjustment  factor  (Ca)  is  used  to  increase  the  cost  of  the  original 
construction.  This  is  done  to  allow  for  an  increase  in  contracting  cost 
to  have  a  site  quickly  repaired. 

Time  of  overtopping  is  determined  as  follows:  Consider  a  bridge 
site  located  at  a  stream  which  has  a  stage-discharge  curve,  I,  as  shown  in 
Figure  7.  A  bridge  having  an  opening  size  of  L]  is  constructed  over  the 
stream,  resulting  in  a  new  stage-discharge  curve.  Notice  that  curve  II  is 
obtained  using  the  assumption  that  no  overtopping  occurs,  i.e.,  H  =  ». 
Now  assume  a  design  scheme  which  has  a  combination  of  bridge  length  L]  and 
an  embankment  height  of  H-| .  Then  Qe  is  the  peak  discharge  producing  a 
water  surface  elevation  equal  to  the  height  of  the  embankment  (point  B). 
Any  flood  producing  a  peak  flow  larger  than  Qe  will  overtop  the  embankment. 

The  time  of  overtopping  for  each  peak  flow  is  determined  from 
the  individual  hydrograph,  e.g.  t2>  t3,  t4  in  Figure  7.  These  values  are 
then  used  to  compute  embankment  erosion  and  traffic-related  losses. 

3.     Bridge  Pier  Scour 

In  major  floods,  bridges  frequently  fail  because  of  scour  around 
piers  and  abutments.  Depth  of  scour  is  affected  by  both  hydraulic  and 
geologic  conditions,  including  type  of  substrata,  pier  shape,  angle  of  the 
pier  to  the  flow,  protective  measures,  depth  of  flow,  velocity  and  turbulent 
intensity.  Among  the  more  random  variables  is  debris  clogging,  which  has  been 
known  to  increase  scour  depth  significantly  and  cause  unexpected  failure. 

36 


4> 

I 

c 

C 

c 
a> 

0) 

e 

8 

^ 

at. 

^^> 

c 

a> 

o 

o> 

c 

JO 

■o 

9 

E 
Ll) 
ii 

ds 

it 

£  x 

-I 

UJ 

o 

z 

> 

' 

E 


CD 
C 
•r" 

Q. 

Q- 
O 
+-> 

S- 

> 
o 

+-> 

E 


CO 

-Q 

E 

LU 

M- 
O 

E 

o 

•I- 

+J 

e 


s- 

<L> 
Q 


37 


Methods  for  computing  scour  depth  are  still  evolving.  The  existing 
formulas  are  derived  primarily  from  laboratory  experiments  and  while  these 
formulas  predict  scour  depth  qualitatively,  they  do  not  incorporate  the 
salient  time-varying  character  of  scour.  Since  the  laws  governing  sediment 
transport  have  not  been  adequately  defined,  model  results  using  presently 
available  equations  may  not  reproduce  the  prototype  behavior.  In  an  approach 
to  this  problem,  the  FHWA  is  sponsoring  a  massive  field  program  to  collect 
scour-related  data. 

For  modeling  purposes  of  this  study,  the  empirical  equation  by 
Laursen  has  been  adapted  for  scour  computation.  The  model,  however,  is 
sufficiently  flexible  to  accept  a  more  realistic  scour  formula  should  one 
become  available.  The  following  simplified  assumptions  are  implemented 
in  the  model  for  scour  analysis  without  compromising  accuracy  to  a  signifi- 
cant degree: 

1.  Pier  scour  is  the  same  for  all  shapes  of  piers,  and 

2.  Velocity  across  the  bridge  opening  is  uniform. 

The  scour  analysis  implemented  allows  the  user  three  different  options: 

1.  Footing  is  set  on  solid  rock.  This  implies  no  scour 
loss. 

2.  Piles  rest  on  solid  rock.  This  is  taken  as  a  cost 
added  to  the  initial  construction  cost. 

3.  Piers  rest  on  alluvial  soil. 

The  first  and  second  options  imply  that  the  scour  risk  due  to  flooding 
is  nil.  The  implementation  of  the  third  option  makes  use  of  the  follow- 
ing empirical  depth-of-scour  equation  presented  by  Laursen  (15). 
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bp      ds 
—  =  5.5  — 

y     y 


i   ds   \  1-7 
+  i    -i 


,11.5  y 


(16) 


VJ 


Where       J)p  =  width  of  rectangular  pier, 
y  =  depth  of  flow,  and 
ds  =  depth  of  scour. 

This  relationship  is  shown  in  Figure  8. 

If  the  computed  depth  of  scour  exceeds  the  user  supplied  maximum 
depth  for  option  3,  bridge  failure  occurs.  In  this  case  the  economic  loss 
is  twice  the  cost  of  the  bridge  times  the  flood  recurrence  interval  since 
the  bridge  is  lost  and  a  new  one  must  be  built  in  its  place. 

Scour  which  does  not  cause  complete  bridge  failure  is  assumed  to 
cause  some  increase  in  losses  determined  by  multiplying  the  scour  depth 
by  a  per  unit  cost  to  repair  the  scour  damage.  Considering  that  there  is 
little  documentation  by  the  dynamic  formation  or  effects  of  nondestructive 
scour,  this  method  provides  a  reasonable  estimate  for  use  by.design_  ,:, 
engineers. 

TRAFFIC-RELATED  LOSSES 

A  major  portion  of  the  economic  loss  resulting  from  damage  to  a 
stream  crossing  structure  is  the  traffic  stoppage  or  the  delay  caused  by 
an  inconvenient  detour.  The  total  time  that  the  traffic  is  not  allowed  to 
travel  at  its  normal  rate  over  the  crossing  is  assumed  to  be  equal  to  the 
sum  of  the  flood  overtopping  duration  and  the  additional  time  required  to 
repair  significant  damage  to  the  site.  The  duration  of  overtopping  is 
computed  in  the  RISK  module.  The  time  of  repair  is  estimated  from  a  graph 
similar  to  Figure  9  which  is  developed  by  the  analyst  for  each  case. 
The  distribution  and  magnitude  of  the  average  daily  traffic  across  the 
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Figure  9.  Assumed  Traffic  Restoration  as  a  Function 
of  Embankment  Erosion  (13) 
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bridge  is  also  required  for  each  case.  There  are  four  sub-categories  of 
traffic-related  losses: 

1.  Increased  running  cost  due  to  detour, 

2.  Lost  time  of  vehicle  occupants  due  to  increased  time  of  travel, 

3.  Accidents  on  additional  detour  miles,  and 

4.  Accidents  due  to  the  unexpected  obstacle. 

After  barricades  are  placed  around  the  site,  the  first  three  types  of  losses 
are  incurred  on  the  detour.  The  fourth  type  is  the  cost  of  the  accidents 
expected  at  the  stream  crossing  and  is  postulated  to  occur  as  traffic  comes 
upon  the  failure  immediately  after  the  failure  occurs. 

Increased  running  cost   is  the  difference  between  running  cost  on 
the  detour  and  the  normal  route.  This  cost  is  computed  as  a  function  of 
average  daily  traffic,  travel  distance,  duration  of  detour,  design  speed, 
and  vehicle  distribution.  Detour  duration  is  the  sum  of  the  overtopping 
duration  plusythe  repair  time.  Overtopping  duration  is  computed  in  the 
flood  routing  calculations.  Five  basic  classes  of  vehicles  are  assumed 
and  actual  traffic  distributions  are  fit,  as  closely  as  possible,  to  the 
five  basic  classes  (13).  These  are: 

1.  0.7-ton  passenger  cars, 

2.  1.25-ton  commercial  delivery  vans, 

3.  1.55-ton  single  unit  trucks, 

4.  2.2-ton  gasoline  semitrailer  trucks,  and 

5.  2.75-ton  diesel  semitrailer  trucks. 

Running  costs  are  estimated  for  passenger  cars  on  zero  grades  and  these 
costs  are  adjusted  to  reflect  the  other  classes  of  vehicles  in  the  five- 
class  distribution.  The  passenger  car  running  costs  as  a  function  of  speed 
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are  fit  to  a  parabola  as  a  means  of  facilitating  computer  applications.  The 
low  point  of  the  parabola,  or  most  economical  running  speed,  occurs  between 
30  and  40  miles  per  hour. 

Time  losses   are  a  function  of  average  daily  traffic,  detour  duration; 
travel  distance,  vehicle  occupancy  rate,  design  speed,  and  the  value 
placed  on  an  individual's  time.  The  occupancy  rate  and  time  values  are 
averages  which  apply  to  all  individuals  in  the  five  classes  of  vehicles. 
Only  the  difference  in  time  between  the  normal  route  and  the  detour  is 
considered. 

The  death  rate  is  used  as  the  basic  unit  to  compute  the  increased 
aeeident  costs   imposed  by  the  detour  (13).  For  each  death,  there  are  a 
certain  number  of  personal  injuries  and  a  certain  number  of ^property  damage 
accidents,  each  of  which  may  or  may  not  be  associated  with  a  death.  For 
example,  30  personal  injuries  for  each  traffic  death  and  300  property 
damage  accidents  might  occur  for  each  death.  The  accident  losses   are 
computed  by  applying  the  above  rates  to  costs  of  deaths,  personal  injuries 
and  property  damage.  These  losses  are  computed  on  a  vehicle  mile  basis-- 
the  standard  for  death  rate  statistics.  Accident  losses  are  computed  as 
a  function  of  average  daily  traffic,  length  of  detour,  duration  of  detour, 
death  rate,  ratio  of  personal  injuries  to  deaths,  cost  of  death,  cost 
of  personal  injury,  and  cost  of  a  property  damage  accident. 

The  last  loss  category  is  the  expected  accident  due  to  the 
unexpected  obstacle   at  the  stream  crossing  site.  The  higher  death  rate  is 
defined  as  the  product  of  the  death  rate  for  normal  conditions  over  one 
mile  of  roadway,  and  a  death  rate  factor  for  unexpected  obstacles  (13). 
This  death  rate  factor  is  somewhat  subjective  because  of  the  lack  of  data 
on  this  type  of  accident.  Perhaps  a  factor  of  one  thousand  is  appropriate. 
The  engineer  can  also  modify  the  ratios  of  personal  injuries  and  property 
damage  accidents  to  deaths  since  these  ratios  may  differ  from  those  on 
the  detour. 
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The  parameters  necessary  to  evaluate  these  losses  are: 

X]  =  Duration  of  Detour  =  Duration  of  Overtopping  +  Repair  Time  (hrs) 

X2  =  Average  Daily  Traffic,  ADT  (Vehicles/Day) 

X3  =  Passenger  Cars  (Fraction  of  ADT) 

X4  =  Commercial  Delivery  Vehicles  (Fraction  of  ADT) 

X5  =  Single  Unit  Trucks  (Fraction  of  ADT) 

X5  =  Gasoline  Service  Trailers  (Fraction  of  ADT) 

X7  =  Diesel  Semi -Trailers  (Fraction  of  ADT) 

X8  =  Length  of  Detour  (Miles) 

Xg.i  =  Normal  Distance  (Miles) 
Xg  =  Speed  on  Detour  (Miles/Hr) 

Xgj  =  Normal  Speed  (Miles/Hr) 
Xio  =  Occupancy  Rate  (People/Vehicle) 

Xn  =  Accident  Distribution  Ratio  -  Normal  Conditions  (Personal 
^In  juries/Death) 

X]2  =  Accident  Distribution  Ratio  -  Normal  Conditions  (Property 
Damage/Death) 

X13  =  Accident  Distribution  Ratio  -  Unexpected  Obstacle  (Personal 
Injuries/Death) 

X-|4  =  Accident  Distribution  Ratio  -  Unexpected  Obstacle  (Property 
Damage/Death) 

X-J5  =  Death  Rate  (People/100  Million  Miles) 

X-J6  =  Death  Rate  Factor  for  Unexpected  Obstacle  (Multiplier  to  X-^) 

C-|  =  Cost  of  Death  ($) 

C2  =  Cost  of  Personal  Injury  ($) 

C3  =  Cost  of  Property  Damage  ($) 

C4  =  Value  of  Time  ($/Hr) 

C5  =  Running  Cost  of  a  Passenger  Car  over  the  Detour  Distance 
($/100  Vehicles) 
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C5  =  Running  Cost  of  a  Passenger  Car  over  the  Normal  Distance 
($/100  Vehicles) 


Parameters  X]  through  X-|0  are  different  for  each  site  considered,  but  Xn 
through  C3,  except  X-j 5 „  represent  national  statistics.  X]g  is  a  multiplier 
applied  to  the  death  rate  to  allow  for  the  increased  hazard  of  an  unexpected 
obstacle.  The  value  of  the  parameter  varies  depending  on  the  site  conditions, 
The  value  of  time  (C4)  is  based  on  the  average  income  per  year  of  the 
users  of  the  bridge.  The  results  tabulated  below  (Table  2)  are  developed 
from  a  study  performed  by  the  Stanford  Research  Institute  (16).  The 
results  shown  in  the  table  are  based  on  the  report-adjusted  data  set  with 
a  sample  size  of  807  using  the  equation 


C4  =  1.803  +  0.4611 


(17) 


where 


I  =  income  level . 


Table  2.  Value  of  Time  ($/hr)  by  Income  of  Highway  User 


I  =  1 


I  =  2 


I  =  3 


I  =  4 


I  =  5 


I  =  6 


I  =  7 


I  <=   8 


C4   2.26 
where 


2.73 

=  1 

=  2 

=  3 

=  4 

=  5 

-  6 

=  7 


3.19 


3.64 


4.11 


4.57 


5.03 


s  income  level  under  $4,000/year, 

s  income  between  $4,000  -  5,999/year, 

s  income  between  $6,000  -  7,999/year, 

s  income  between  $8,000  -  9,999/year, 

s  income  between  $10,000  -  11,999/year, 

s  income  between  $12,000  -  14,999/year, 

s  income  between  $15,000  -  20,000/year,  and 


=8  is  income  over  $20,000/year. 


5.49 
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The  running  cost  of  a  passenger  car  in  dollars  per  1000  vehicle 
miles  is  a  function  of  speed  as  shown  in  Figure  10.  The  appropriate  . 
equations  for  C5"  and  C5  are 

C5  =  42.5  -  .455X9  +  -0068X92  (18) 

C6  =  42.5  -  .455X9 j  +  .0068X9>12 

To  adjust  these  passenger  car  running  costs  for  varying  types  of  vehicle 
distributions,  Equation  18  becomes 

C5  =  (42.5  -  .455Xg  +  .0068Xo2  (19) 

•  (X3  +  1 .2X4  +  2.OX5  +  3-2X6  +  3JX7 
and               C6  =  (42.5  -  .455Xg  +  .0068Xo2) 

•  (X3  +  1.2X4  +  2.OX5  +  3-2x6  +  3.IX7) 

To  estimate  the  losses  associated  with  running  costs,  it  is 
necessary  to  compute  the  running  cost  over  the  normal  route  and  the  detour. 
The  difference  in  these  two  cost  values  represents  the  additional  cost  to 
the  user  of  having  to  detour  due  to  the  failure  of  the  bridge.  The  equation 
for  computing  the  running  cost  ($)  is: 

xl  x2  x8  c5    xl  x2  x8.1  c6  ,  t 

XL]  = (20) 

24,000        24,000 

The  time  lost  by  the  vehicle  occupants  is  the  delay  or  the  additional  time 
it  takes  to  detour  the  site.  The  value  of  lost  time  ($)  is  computed  by 
calculating  the  difference  of  the  time  value  of  the  detour  and  the  original 
route.  The  equation  used  for  this  calculation  is: 

8    8  1   2 
XL2  "  xl  x£-  X^7  24  X10  C4  (21) 
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Figure  10.  Passenger  Car  Running  Costs  (17) 
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The  expected  accident  cost  ($)  due  to  the  difference  in  the  dollar  yalue 
of  accidents  on  the  detour  and  on  the  original  route  is  calculated  by: 

Xl  X0      (Xp  -  Xp  1  )  X-i  r 

XL3  =  -1— ! * 5 |d !£     (C     +  xnCo  +  X1?CJ  (22) 

J       2.4  x  109  '    "  L  U   J 

The  expected  accident  cost  due  to  an  unexpected  obstacle  is  computed 
by  assuming  one  mile  of  road  has  a  one-hour  exposure  to  a  higher  death 
rate  defined  as  the  death  rate,  X15,  times  a  death  rate  multiplier  for 
unexpected  obstacles,  X-jg.  Thus,  the  higher  rate  is  X-|5X,,-.  The  accident 
distribution  ratios  are  X13  and  X-j 4  which  may  vary  from  those  for  the 
normal  death  rate,  X-|-j  and  X-| 2.  The  equation  for  calculating  this  loss 
($)  is: 

XL4  =  2  15  I6     (Cl  +  X13C2  +  X14C3)         (23) 
4   2.4  x  10y 

By  applying  Equations  18  through  23  to  a  set  of  data  describing 
a  given  site,  it  is  possible  to  compute  the  total  dollar  value  of  the 
traffic  loss  due  to  flooding  at  the  culvert  site.  Table  3  shows  the 
traffic-related  data  for  two  types  of  highway  systems. 

FLOOD  LOSSES 

The  constriction  of  waterways  due  to  approach  embankments  incurs 
incremental  damage  to  the  adjacent  properties  when  severe  flooding  occurs. 
This  incremental  damage  is  an  amount  of  loss  in  addition  to  the  damage 
that  would  be  caused  before  the  embankments  are  constructed.  The  magni- 
tude of  the  damage  depends  on  the  specific  land  uses  and  the  type  of 
development  occurring  in  the  flood  plains,  for  example,  woodlands,  pasture 
lands,  farm  crops,  private  homes,  industries  or  any  combination  of  these. 
Different  levels  of  damage  are  associated  with  these  land  uses  by  flooding. 
For  example,  pasture  and  woodlands  may  be  subject  to  little  damage  from 
water,  while  a  private  home  or  industry  may  be  destroyed  by  severe  flooding. 
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Table  3.     Traffic  Loss  Data 


Variable'  Four-Lane  Rural  Units 

Highway  Highway 


h 

16000 

h 

85 

h 

.01 

h 

.02 

h 

.03 

h 

.09 

X8"X8.1 

1.20 

X9 

55 

X10 

1.7 

Xn8 

30 

v   * 

X12 

300 

V   * 

X13 

15 

X14 

150 

Y   * 

X15 

5.5 

lo 

1000 

r   * 
Ll 

50000 

c2* 

2000 

L3 

400 

L4 

2 

*Accident  Facts, 

Nation 

466 

Vehicles/day 

.955 

Fraction  of  ADT 

.043 

Fraction  of  ADT 

0 

Fraction  of  ADT 

.002 

Fraction  of  ADT 

0 

Fraction  of  ADT 

1.41 

Mi  1  es 

25 

Miles/hr 

2 

People/Vehicle 

30 

Personal  Injury 
Incidences/Death 

300 

Property  Damage 
Claims/Death 

15 

Personal  Injury 
Incidences/Death 

150 

Property  Damage 
Claims/Death 

5.5 

People/ 100  Million 

miles 

500 

50000 

$ 

2000 

$ 

400 

$ 

2 

$/hr 
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Flood  loss  is  usually  expressed  by  a  stage-damage  curve  for  the 
area  inundated  by  the  flood.  A  technique  for  evaluating  flood  damage  at 
culvert  sites  was  developed  by  WRE  (13).  This  method  relates  the  flood 
loss  to  the  ponding  of  flood  water  upstream  from  the  culvert  site.  The 
computational  procedure  is  similar  to  the  one  generally  used  in  reservoir 
projects.  In  either  case  the  assumption  of  horizontal  pool  level  has  been 
made.  Thus  the  stage-area,  and  hence  the  stage-damage,  curve  can  easily 
be  developed. 

At  a  bridge  site  the  situation  is  somewhat  different.  Instead  of 
having  a  horizontal  water  level,  the  water  surface  elevation  varies  with 
flood  magnitudes  as  well  as  with  location.  To  develop  a  single  stage- 
damage  curve  under  such  conditions  is  virtually  impossible. 

Theoretically  a  bridge  owner  should  not  be  held  responsible  for 
flood  losses  incurred  under  natural  flow  conditions  before  a  bridge  is 
in  place,  no  matter  how  large  the  flood  is.  It  is  the  bridge  owner's 
liability,  however,  for  the  additional  portion  of  flood  damage  which  is 
caused  by  the  backwater  associated  with  the  construction  of  the  stream 
crossing  structures.  The  stage-damage  function  developed  in  this 
study  is  based  upon  this  concept. 

The  first  step  in  the  process  of  deriving  the  stage-damage  curves  is 
defining  the  steady  state  hydraulic  regime.  As  discussed  previously, 
water  surface  elevations  in  the  vicinity  of  the  site  (assuming  no  bridge 
is  in  place)  are  calculated  by  the  two-dimensional  computerized  program 
FEM.  Figure  11a  shows  the  water  surface  profiles  computed  for  natural 
conditions  for  flows  Q-, ,  Q~  and  Q~.  The  next  step  is  to  calculate .. 
water  surface  elevations  assuming  the  bridge  is  in  place  (Figure  lib). 
The  constriction  will  cause  an  increase  in  water  surface  elevation  in 
the  vicinity  of  the  bridge  crossing.  Increased  water  surface  elevations 


50 


a)     W.S.   PROFILES    WITH    NO  BRIDGE 

BRIDGE 


b)     W.S.   PROFILES    WITH    BRIDGE   IN    PLACE 


LOSS    ($) 
C )     STAGE  -  DAMAGE    DIAGRAM 


Figure  11.     Stage-Damage  Curves 
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result  in  increased  flood-related  damages  due  to  the  areal  spreading  of 
the  flood  and  additional  damages  at  locations  affected  with  no  bridge 
in  place.  The  two  sets  of  surface  elevations  is  the  starting  point  for 
estimating  the  economic  losses  associated  with  a  flood  of  known  magnitude 
both  for  the  natural  conditions  and  the  altered  circumstance  of  having 
the  bridge  in  place.  The  difference  of  these  two  economic  losses  defines 
the  economic  loss  associated  with  a  given  bridge  design  for  a  specific 
flood. 

For  a  specific  site  the  two  hydraulic  regimes  (with  and  without 
the  bridge)  are  plotted  on  a  USGS  contour  map  of  the  area.  This  defines 
the  areal  coverage  and  actual  water  depth  for  both  regimes. 

Once  the  hydraulic  regimes  are  defined  for  the  cases  with  and 
without  the  bridge,  the  calculation  of  the  economic  losses  becomes 
straightforward.  Flood  damage  tables  similar  to  Table  4  are  generated 
for  any  field  site  that  is  to  be  analyzed.  The  details  of  the 
development  of  these  tables  can  be  found  in  Appendix  B  of  Reference  13. 

With  the  required  tables,  the  hydraulic  information  and  a  knowledge 
of  the  land  use  pattern  in  the  vicinity  of  the  proposed  bridge  opening, 
it  is  possible  to  calculate  the  flood-related  losses  for  the  site  with 
and  without  the  bridge.  The  difference  in  these  values  is  the  flood 
damage  due  to  the  constriction  caused  by  the  bridge  crossing.  Figure  lie 
shows  stage-damage  curves  derived  from  Figure  lib. 

In  this  study  these  calculations  for  the  economic  damage  caused 
by  the  backwater  are  computed  manually  using  the  computer-derived 
hydraulic  information.  Data  are  computed  for  the  complete  array  of  decision 
parameters  considered,  various  discharges,  embankment  heights  and  bridge 
lengths.  The  resulting  values  are  part  of  the  input  data  read  into  the 
risk  program  model  and  used  in  determining  the  total  economic  loss  and 
resulting  risks  for  all  cases  of  interest. 
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RISK  ANALYSIS 

Risk  analysis  is  defined  as  the  evaluation  of  the  losses 
incurred  while  playing  strategy  i  due  to  a  series  of  possible  states 
of  nature  and  an  assessment  of  the  probability  distribution  of  the 
states  of  nature. 

We  define  the  risk  for  strategy  i,  R-j ,  for  a  discrete  proba- 
bility distribution  of  the  states  of  nature  as: 

n 

*i  =E  Yj  <24> 

j=i 

n 
with        Y]     P.  =  1 .0 
j=l 


where       P.  =  probability  of  occurrence  of  state  of  nature  j,  and 
j 

D-  =  loss  incurred  by  the  player  due  to  the  occurrence  of 


the  j   state  of  nature, 


For  a  continuous  distribution  of  the  state  of  nature,  we  define 
the  risk  R.  for  strategy  i  as: 


oo 


R,  =  JL(Q)f(Q)dQ  (25) 

0 


with        /f(Q)dQ  =  1,  f(Q)  >  0 
0  0<Q<c° 
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where       Q  i- 


f(Q)   defines  the  probability  density  function  of 
the  states  of  nature  Q,  and 

L(Q)  is  the  continuous  loss  function. 

In  the  case  of  applying  Risk  Analysis  to  select  design  schemes, 
i.e.,  the  combinations  of  bridge  length  and  height  of  approach  embankment, 
we  make  the  following  analogies  to  the  theory  expressed  above.  The 
player's  strategy  i  becomes  (design  scheme).,  which  in  turn  is  determined 
by  two  design  parameters:  (embankment  height).,  and  (bridge  length).. 
These  two  design  parameters  are  referred  to  as  decision  variable  set  i. 
The  player  corresponds  to  society  in  general.  The  possible  states  of 
nature  correspond  to  the  possible  floods  at  a  bridge  site.  We  note  that 
the  distribution  of  possible  floods  at  the  bridge  site  is  a  continuous 
distribution  and  can  be  approximated  by  a  double  exponential  or  Gumbel 
distribution.  Therefore  we  must  consider  Equation  25  as  the  governing 
risk  equatioji  which  must  be  evaluated  for  each  design  scheme  i,  to 
determine  that  scheme's  total  economic  response. 

The  RISK  Model  approximates  the  integral  of  Equation  12  by 
using  the  five  flood  events,  Q-, ,  Q2,  Q3»  Q4  and  Q5  as  shown  below: 

Qi       Q2 

Ri  =yL(Q)f(Q)dQ  =yL(Q)f(Q)dQ  +yL(Q)f(Q)dQ 
0  0  Ql 

Q3  Q4         ^5 

+  /L(Q)f(Q)dQ  =/L(Q)f(Q)dQ  +/L(Q)f(Q)dQ 
Q2         Q3         Q4 
00 
+  y[(Q)f(Q)dQ  (26) 

Q5 


56 


Each  of  the  six  terms  on  the  right  hand  side  of  Equation  26  is  then 
approximated  by  the  following  expressions: 

Ql 

jL(Q)f(Q)dQ  »  0    since  it  is  assumed  L(Q)  =  0, 
0  for  0<Q<Q1. 


(27a) 


Q2  Q2 

y*L(Q)f(Q)dQ   «  L2/f(Q)dQ  =  L2|p[Qi<Q<Q2 
Ql  Ql 


(27b) 


Q3  Q3 

/L(Q)f(Q)dQ  «  l3ff(Q)dQ   =  L3  p[q2<Q<Q3 
Q2  Q2  ( 


(27c) 


Q4  Q4 

JL(Q)f(Q)dQ  «  LAff(Q)dQ   =  L4  P 

Q3        Q3 


Q3lQ<Q4 


(27d) 


Q5 

/L(Q)f(Q)dQ 
Q4 


Q5 
L5/(Q)dQ 
Q4 


L5  P|Q4<Q<Q5 


(27e) 


00  00  , 

JL(Q)f(Q)dQ  «  L5yf(Q)dQ  =  L5  P  [q>Q5 

is  "5 

since  it  is  assumed  L(Q)  =  1$,  for  Q>Q5 


(27f) 


where 


Li ,  L2,  L3,  L4  and  L5  are  the  losses  associated  with  Q] , 
Q2>  Q3S  Q4  and  Q5,  respectively  for  a  given  design  scheme 
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Probability,  P,  in  Equations  27a  through  27f,  weights  the  flood  loss  in  each 
flood  interval.  The  application  of  this  concept  for  risk  analysis  will 
be  further  illustrated  in  Chapter  V,  with  specific  examples. 


SELECTION  OF  MOST  ECONOMIC  DESIGN 

The  sum  of  the  estimated  construction  costs  of  a  design  scheme 
and  its  calculated  flood  risks  equals  the  economic  response  of  that 
design  scheme.  The  design  goal  of  a  highway  stream  crossing  is  a  design 
which  minimizes  the  economic  response.  This  process  involves  a  great 
deal  of  computational  effort  by  using  varying  decision  variables,  one 
at  a  time.  The  computation  may  become  more  complicated  as  the  number 
of  decision  variables  increases. 

In  the  design  of  highway  stream  crossings  the  two  most  important 
decision  variables  are  the  size  of  the  bridge  opening  and  the  height 
of  the  approach  embankments.  Other  decision  variables  related  to  scour 
around  piers  and  abutments,  though  important,  are  not  included  in  the 
analysis  procedure  described  in  this  report  because  these  parameters  have 
not  yet  been  adequately  defined.  These  parameters  should  be  incorpor- 
ated in  the  present  model  when  better  information  becomes  available. 
Decision  variables  regarding  protection  measures  on  embankments  are 
also  omitted  in  the  analysis.  Consideration  should  also  be  given  to 
include  these  variables  into  the  model  in  future  studies. 

For  illustration  purposes,  assume  that  there  are  an  infinite 
number  of  design  schemes  consisting  of  an  infinite  number  of  combinations 
of  the  bridge  decision  variables,  embankment  height  (H)  and  bridge 
length  (L).  Further,  assume  L,H  >  0.  For  each  such  combination  of 
the  decision  variables  L  and  H  we  evaluate  the  risk  according  to 
Equation  25,  as  f(L,H).  We  next  determine  the  construction  cost  for  the 
design  scheme  corresponding  to  decision  variables  embankment  height  H 
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and  bridge  length  (L)  as  C(L,H).  The  total  economic  response  function, 

TER(L,H),  is  defined  in  terms  of  the  decision  variable  embankment  height 

(H),  and  bridge  length  (L)  as: 

Construe- 
Risk    tion  cost 

TER(L,H)  =  f(L,H)  +  C(L,H)  ,  L,H  >  0  (28) 

This  response  surface  is  sketched  in  Figure  12.  The  optimal  design 

corresponds  to  that  choice  of  the  decision  variables  L  .  and  H  .  which 

opt     opt 

correspond  to  the  minimum  value  of  TER(L,H);  namely  I  =  TER  (L  .,  H  .), 

opt   opt  • 

as  shown  in  Figure  12. 

This  study  considers  a  finite  number  of  design  schemes  involving 
three  embankment  heights  and  three  bridge  lengths.  Then  the  model  evaluates 
the  total  economic  response  for  each  of  the  seven  possible  design  combi- 
nations. 
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IV.  APPLICATION  OF  RISK  MODEL  TO  AN  EXAMPLE  PROBLEM 


SITE  INFORMATION 

The  principles  of  analysis  presented  in  Chapter  III  are  applied 
to  the  Tall  anal  la  Creek  bridge  site  near  Waldrup,  Mississippi.  The 
flood  plain  consists  of  woodlands  and  parks.  At  the  bridge  site  the 
drainage  basin  of  Tallahalla  Creek  is  approximately  33  square  miles. 
The  site  map  is  shown  in  Figure  13.  At  present  there  is  a  500-foot 
bridge  across  Tallahalla  Creek  near  Waldrup.  The  details  are  shown  in 
Figure  14. 


DESIGN  ALTERNATIVES 

In  the  analysis  seven  alternative  bridge  design  schemes  are 
considered  including  the  present  bridge.  These  design  schemes  are 
listed  in  Table  5. 

Table  5  shows  that  the  two  major  decision  variables  considered 
are  the  bridge  length  (or  bridge  opening  size)  and  the  height  of  approach 
embankments.  While  these  two  variables  are  the  primary  factors  govern- 
ing the  overall  economic  picture  at  the  crossing  site,  other  factors, 
such  as  type  of  stream  crossing  (e.g.  skewed,  eccentric,  normal,  or 
any  of  the  combinations),  bridge  type  (e.g.  slab,  I-beam  composite,  or 
steel  bridge),  and  the  class  of  highway  (e.g.  primary  or  secondary), 
also  affect  the  result  of  economic  analysis.  The  type  of  crossing 
governs  the  amount  of  backwater  and  thus  property  damage  and  traffic 
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SCALE 


Figure  13.  Site  Map  at  Tallahalla  Creek,  Mississippi 
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delay.  Bridge  type  determines  the  unit  cost  of  the  bridge,  while  the  class 
of  highway  affects  the  specification  of  design  standards  for  any  given 
bridge  length  and  embankment  height.  The  extent  to  which  these  additional 
factors  would  influence  the  economic  analysis  depends  on  each  individual 
bridge  site.  For  the  present  analysis  it  is  assumed  that  these  additional 
factors  be  held  constant.  That  leaves  the  economic  analysis  to  be  focused 
upon  the  two  degrees  of  freedom,  i.e.,  the  bridge  length  and  the  embankment 
height. 


Table  5.   List  of  Design  Alternatives 


Decision  Variables 

Design  Alternative  Bridge  Length  (ft)    Embankment  Elevation  (ft) 

1  300  315 

2*  500  315 

3  500  310 

4  500  320 

5  1000  315 

6  1000  310 

7  1000  320 


• 


This  scheme  is  the  present  design  in  place 


HYDRAULIC  COMPUTATIONS 

Hydraulic  computations  for  all  the  design  alternatives  were  per- 
formed using  the  Finite  Element  Model  presented  in  Volume  III  of  the  final 
report  of  this  project.  The  network  system  used  to  represent  the  flow 
region  for  the  500-foot  bridge  opening  is  the  same  as  that  used  in  the 


64 


example  cases  of  the  Tallahalla  Creek  site  in  Volume  III.  Minor  modification 
to  the  network  system  across  the  bridge  openings  of  300  and  1000  feet 
was  required  to  reflect  the  change  in  opening  size.  Pertinent  data  for 
the  opening  configuration  are  shown  in  Figure  15  and  Tables  6  through  8.. 


Table  6.   Embankment  Data  for  Bridge  1 
(Bridge  Length  =  300  ft) 


Node 


Ground  Elevation 


Length 


Top  Elevation 


158 
147 
136 
125 
114 
104 
115 
126 


309.0 

0.0 

315.0 

307.7 

365.0 

315.0 

306.5 

365.0 

315.0 

305.7 

220.0 

315.0 

305.0 

220.0 

315.0 

300.0 

320.0 

315.0 

304.0 

330.0 

315.0 

308.0 

330.0 

315.0 
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Table  7.   Embankment  Data  for  Bridge  2 
(Bridge  Length  =  500  ft) 


Node      Ground  Elevation       Length       Top  Elevation 


158 
147 
136 
125 
114 
104 
115 
126 


309.0 

0-0 

315.0 

307.7 

365.0 

315.0 

306.5 

365.0 

315.0 

305.5 

120.0 

315.0 

304.6 

120.0 

315.0 

300.0 

520.0 

315.0 

304.0 

330.0 

315.0 

308.0 

330.0 

315.0 

158 
147 
136 
125 
114 
104 
115 
126 


Table  8.   Embankment  Data  for  Bridge  3 
(Bridge  Length  =  1 ,000  ft) 


Node      Ground  Elevation       Length       Top  Elevation 


309.0 

0.0 

315.0 

307.5 

265.0 

315.0 

306.0 

265.0 

315.0 

306.3 

107.5 

315.0 

306.5 

107.5 

315.0 

307.0 

1020.0 

315.0 

307.5 

192.5 

315.0 

308.0 

192.5 

315.0 
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HYDROLOGY 

The  flood  frequency  curve  for  the  study  site  is  shown  in  Figure  4 
(page  26).  For  each  bridge  design,  it  is  necessary  to  consider  a  wide  range 
of  hydrologic  possibilities.  Since  flow,  Q,  is  considered  a  continuous 
random  variable,  it  is  improper  to  speak  of  the  probability  associated 
with  one  particul'ar  flood  event.  Instead,  one  must  consider  the  probability 
that  a  storm  event  will  lie  in  the  closed  interval  [a,  b]  as  shown  in  the 
equation  below. 

b 
P  [a<Q<b]  =  /  f(Q)  dQ  (29) 

a 


where 


f(Q)  =  the  probability  density  function  of  the  random 
variable,  Q. 


By  referring  to  Figure  4  one  may  determine  the  recurrence  interval, 
r. ,  for  a  givefi  flow,  Q  .  The  relationship  between  r.  and  Qr  is 
expressed: 


where 


P  (Q>Qr)  =  p.  (30) 

Q  =  flow, 

Q  =  flow  associated  with  recurrence  interval,  r . ,  and 
r.  =  recurrence  interval  associated  with  flow,  Q  . 

Thus  the  flow  associated  with  a  storm  having  a  recurrence  interval  of 
ten  years  (e.g.,  a  ten-year  storm)  is  exceeded  10  percent  of  the  years. 


68 


For  the  analysis  performed  in  the  RISK  module,  five  storm  events 
are  considered  as  shown  in  Table  9.  It  was  assumed  that  storms  with  a 
peak  discharge  below  15,000  cfs  caused  zero  weighted  risks  and  that  those 
above  27,500  cfs  did  not  cause  weighted  risks  above  that  caused  by  the 
200-year  storm  event.  The  situation  is  shown  in  Figure  16.  Flow  Q  is 
a  random  variable  with  a  probability  density  function,  f(Q). 


Table  9.  Storm  Events 


Probability  that 

Recurrence 

Flow 

P 

h 

<Q<Qn 

R(i+1) 

■*■ 

T3  +-> 

flow  wil  be 

Interval 

(cfs) 

Floo 
Even 

eqi 
in 

aled  or  exceeded 
any  one  year 

(yrs) 

i  =  1...5 

i  = 

1....4 

qri 

.05 

20 

15000 

.01 

«R2 

.04 

25 

16200 

.02 

qR3 

.02 

50 

19700 

.01 

QR4 

.01 

100 

23500 

.005 

QR5 

.005 

200 

27500 

Qn  <Q<Qn 


R. 


i+1 


LOSS  FUNCTION  AND  RISK 

A  loss  function  L(Q)  is  associated  with  the  random  variable  Q. 

It  is  assumed  that  L(Q)  =  0  for  all  flows  below  the  20-year  flood,  Q2QS 

and  that  L(Q)  =  L    for  all  flows  above  the  200-year  flood,  Q9nn. 
max  ^uu 
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The  Risk,  R,  is  defined  as: 

Q20  Q200 

00  oo 

R  =  Jf(Q)L(Q)dQ  =  /f(Q)L(Q)dQ  +  Jf(Q)L(Q)dQ  +  Jf(Q)L(Q)dQ   (31) 
0  0  Q2Q         Q20Q 

Rj  Rn  RnI 


The  first  term,  R.,  on  the  right  hand  side  of  Equation  31  above  is  zero 
since  L(Q)  =  0  for  every   Q  contained  in  (0,  Q20)'»  ">-e-»  rt  =  °- 

oo  00 

The  third  term,  Rm,  of  (31)  above,  /f(Q)L(Q)dQ  =  Lmax  /  f(Q)dQ  =  f200Lmax 

Q200  Q200 

oo 

since  for  Q>Q200  ,  L(Q)  =  Lmax  ,  and  Jf(Q)dQ  =  f20Q  . 

Q200 

The  evaluation  of  the  middle  term,  Rjj,  is  the  main  issue.     In  the  RISK 
module,  the  following  approximation  is  made  for  R,.: 

R„   -  P[Qr1<^R2]LQr2  +  P[QR2^R3]L(JR3 

+  P[Qr3^VK4  +  P[QR4^QR5]LQR5  (32) 

From  Table  9:     QR]   =  15,000  cfs,  QR2  =  16,200  cfs,  QR3  =  19,700  cfs, 

QR4  =  23,500,  and  QR5  =  27,500  cfs.     LQR2,  LQR3»LQR4  and  LQR5  are  the 

losses  associated  with  flows  QR2,  QR~,  QRZ,  and  QR5  respectively.     Since 

%z  =  Qonn  =   27,500  cfs,   LQDC  =   L 
^R5       ^200  XR5         max 
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Equation  32  above  is  approximated  as  follows: 

R  =  p[<5r1.<Q.<QR2]lQR2  +  P[QR25Q.<qR3]LQR3 

+    P  [Qr3  <<yfy]  LQR4  +  (P  K^m]  +  f 200  )  Lmax  ( 33) 

Referring  to  Table  9  and  Figure  16: 

P[QR1<Q<QR2]  =   -01  P[QR3^R4]    =   '01 

P[QR2^QR3]=   -02  P^R4^R5]    "   -005 

f200  =   -005 

Equation  33  reduces  to: 

R/  *01L16,200  +  •02L19,700  +  •01L23,500  +  -01L27,500 

In  the  RISK  module  the  losses  for  the  16,200  cfs,  19,700  cfs, 

23,500  cfs,  and  27,500  cfs  storm  events  are  evaluated  and  weighted  as 

shown  in  Equation  33  above  to  determine  the  risk  associated  with  a  particular 
bridge  design. 


LOSS  ESTIMATES 

The  total  loss  resulting  from  each  of  the  above  flood  events 
consists  of  the  following  four  losses: 

1.  Backwater  damages, 

2.  Erosion  damage  to  approach  embankments, 

3.  Scour  damage,  and 

4.  Traffic  delay  and  accident  losses. 
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BACKWATER  DAMAGES 

Backwater  damages  are  due  to  inundation  of  areas  of  the  flood 
plain  which  would  normally  not  be  under  water  were  the  given  bridge 
configuration  not  in  place. 

To  develop  a  method  to  evaluate  the  adjacent  property  loss  due 
to  flooding,  it  is  necessary  to  obtain  information  that  relates  damage  to 
the  type  of  land  development  in  the  upstream  flood  plain.  Information 
is  shown  in  Tables  10  through  15  relating  losses  in  terms  of  1974  dollars 
to  the  following  land  uses: 

1 .  Agricultural , 

2.  Manufacturing, 

3.  Single  Family  Residences, 

4.  Retail  Businesses, 

5.  Selected  Services,  and 

6.  Wholesale  Businesses. 

The  development  of  these  tables  is  described  in  Appendix  B  of  Reference  13 

There  are  several  notable  characteristics  of  these  data.  First, 
flood  losses  do  not  always  increase  for  depths  exceeding  four  feet.  For 
agriculture,  the  unit  flood  loss  is  not  expected  to  significantly  exceed 
the  unit  loss  at  a  depth  of  three  feet  on  the  basis  that  most  crops 
do  not  grow  taller  than  this  height.  A  similar  trend  of  limited  unit 
losses  is  expected  for  single  family  residences  and  selected  services. 
Unit  losses  for  manufacturing,  on  the  other  hand,  can  be  expected  to 
increase  for  depths  greater  than  four  feet  because  stock  piles  of  raw 
materials  and  plant  assets  are  commonly  in  excess  of  heights  of  four 
feet.  Similarly,  increasing  unit  losses  are  expected  for  retail  and 
wholesale  businesses.  For  these  latter  land  uses,  extrapolating  available 
data  to  depths  of  six  or  seven  feet  appears  justifiable. 
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Table  10.  Estimated  Flood  Loss  Data  for 
Agriculture  in  Mississippi1 


Direct  and  Indirect  Damages  in  19742 
In  Dollars  Per  Acre  at  Hater  Depths  of: 

0.51      1 '     _V_  V  + 

Mississippi  $28      $36     $50    $61 

Jasper  County  9       12     15     18 


1U.  S.  Department  of  Commerce,  Bureau  of  the  Census,  1969  United 
States  Census  of  Agriculture,  Vol.  1,  Section  1,  Summary  Data, 
pp.  VII,  5  &  6. 

Ibid.,  Section  2,  County  Data,  pp.  249. 

U.  S.  Department  of  Agriculture,  Soil  Conservation  Service, 
Economics  Guide  for  Watershed  and  Flood  Prevention,  March  1964. 

James,  L.  D.,  A  Time-Dependent  Planning  Process  for  Combining 
Structural  Measures,  Land  Use,  and  Flood  Proofing  to  Minimize 
The  Economic  Cost  of  Floods,  Stanford  University,  August  1964. 

2Monetary  data  is  converted  to  a  1974  base  by  applying  average 
annual  compound  interest  rates  earned  on  long-term  U.  S. 
Government  Securities  to  the  basic  data. 
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Another  characteristic  of  the  flood  data  that  should  be  noted  is 
related  to  geographical  application.  The  flood  loss  data  for  agriculture 
(Table  10)  is  developed  for  Jasper  County  and  for  the  State  of  Mississippi, 
whereas  flood  loss  data  on  single  family  residences  (Table  12)  are 
developed  specifically  for  Jasper  County.  The  flood  loss  data  for  the 
remaining  categories,  namely  manufacturing  (Table  11),  retail  business 
(Table  13),  selected  services  (Table  14),  and  wholesale  business  (Table  15), 
are  applicable  to  the  State  of  Mississippi. 

Tables  10  to  15  are  used  in  conjunction  with  an  inventory  of  the 
land  use  in  the  upstream  flood  plain  to  develop  a  stage  flood  loss 
function.  The  approach  is  to  compute  losses  for  land  area  covered  by 
varying  degrees  of  depth  of  flood  waters.  Because  the  flood  plain  at 
the  bridge  site  consists  of  only  woodlands  and  pasture,  an  approximate 
stage-damage  curve  was  developed  to  determine  the  backwater  damages.  For 
a  more  varied  land  use  pattern  in  the  flood  plain,  the  approach  outlined 
above  is  recommended  to  determine  backwater  losses.  The  approximate 
stage-damage  curve  developed  for  the  study  site  is  shown  in  Figure  17 
and  the  flood  damage  information  for  the  flows  considered  in  the  risk 
analysis  is  shown  in  Table  4. 

EROSION  DAMAGES 

Approach  embankments  may  be  eroded  during  the  flood  events 
considered.  The  amount  of  erosion  is  determined  by  the  time  of  over- 
topping and  the  overtopping  head.  In  order  to  calculate  overtopping 
time,  the  duration  of  the  storm  hydrograph  must  be  determined.  The 
pertinent  hydrograph  data  is  shown  in  Table  16.  A  cost  adjustment  factor 
is  used  to  multiply  the  repair  cost  to  indicate  that  emergency  repair  is 
more  costly  than  normal  repair.  In  addition,  an  inspection  time  of 
two  hours  is  assumed  to  inspect  the  roadway  if  any  overtopping  results 
before  opening  the  roadway  to  traffic  again.  The  erosion  information 
is  shown  in  the  bottom  section  of  Table  16. 
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Figure  17.  Stage-Damage  Curve 
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Table  16.  Hydrologic,  Scour  and  Erosion  Data 


MYQROLUGIC  DATA 
LENGTH  OP  LONGEST  wATtKCUUKSt 
ELEVATION  DIFFERENCE' 


17.000   HUES 
$3,000  #EfcT 


SCOUH  INf-OHrtATIUN  (ALL  LENGTHS  IN  FE6TJ 
CHANNEL  TYPE  IS       1,000 
COST  PER  FUOT  UP  SCOOP  18      «0,000 
MAXlrtUM  SCuUK  DEPTH  IS      -0,000 
VELOCITY  NODE  FOR  SCOUR  IS  »0 

RECTANbULAK  kIDTH  0*  PUR  18     •O.OOO 
NUMBtk  OF  PILES  IS  »0 

PILE  DEPTH  IS      »0,000 


EHOSIQN  INFORMATION 
INCIPIENT  YEL(IFT/SEL;  IS     3,5000 
WEIGHT  (LBS/FT3)'UF  SUK  IS   US, 0000 
CUST'AtfJuSTrtgNT  FACTOR  IS     2,0000 
INSPECTION  TIME(HR)  IS     2,0000 
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SCOUR  DAMAGE 

In  the  example  problem  it  is  assumed  that  the  pieces  of  the 
bridge  rest  on  bedrock  so  that  there  is  no  scour  damage.  The  variables 
which  would  be  used  to  calculate  scour  are  set  to  zero.  The  scour 
information  is  shown  in  the  middle  portion  of  Table  16. 

TRAFFIC  DELAY  LOSSES 

During  the  course  of  the  storm  the  bridge  may  become  impassable 
due  to  overtopping  during  which  embankment  erosion  may  occur.  The  sum 
of  overtopping  time,  inspection  time,  and  embankment  repair  time  consti- 
tutes the  total  delay  time  to  motorists.  A  traffic  distribution  is 
assumed  as  shown  in  Table  17  along  with  a  detour  distance  and  travel 
speed  which  is  compared  to  the  normal  route  distance  and  speed.  The 
results  are  used  to  calculate  the  lost  time  of  the  motorist  which  is 
converted  to  an  economic  loss  by  use  of  the  value  of  time  parameter  shown 
in  the  economic  section  of  Table  17.  Additional  losses  are  accident 
and  obstacle  related  damages.  Accident  rate  and  cost  data  are  also 
shown  in  Table  17. 


TOTAL  ECONOMIC  RESPONSE 

In  order  to  determine  the  total  economic  response  for  each  bridge 
design,  the  total  construction  cost  of  each  bridge  design  must  be  added 
to  the  total  weighted  losses.  The  total  construction  cost  is  calculated 
as  an  equivalent  annual  cost  based  on  a  given  project  life  and  interest 
rate  using  the  capital  recovery  factor.  The  total  construction  cost  of 
the  bridge  is  made  up  of  the  cost  of  the  bridge,  the  cost  of  embankments, 
and  the  roadway  costs.  The  bridge  cost  function  and  the  pertinent 
embankment  and  roadway  cost  information  is  shown  in  Table  18. 
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Table  17.  Economic  and  Traffic  Information 


PROJECT  LIFE  (YKS)  s  100,000 

INTEREST  RATE(PER  CENT)  1         7,000 
AVE,  VALUE  OF  TIM£  (S/HK)  2,730 


AVERAGE  DAILY  TRAFFIC  LEVEL  16000,000 

PASSENGER  CAR  TRAFFIC  LEVEL  ".850 

COMMERCIAL  DELIVERY  VEHICLE  LEVEL  ,010 

SINGLE  UNIT  THUCK  LEVEL  ,030 

GASULlNfc  TRUCK  LEVEL  »030 

DIESEL  SEMI  TRAILER  LEVEL  ,090 

NORMAL  TRAVEL  DISTANCE  «,500 

NURHAL  TRAVEL  SPEtO  55,000 

DtTUUR  [RAVEL  DISTANCE  13,000 

DETUUR  TRAVEL  SPEED  (M/'hR)  40,000 

AVE,  VEH,  OCCUPANCY  RATE*  i»700 


ACCIDENT  RATE  DATA 

DEATH  RATt/100  MILLION  MILES  (NORMAL)  5.500 

PERSONAL  INJURIE-S  (NU./UEAIH)  (NORMAL)  30,000 

PnQPERTY  0AM4GE  (N0,/DEATMJ  (NORMAL)  300,000 

PERSONAL  INJURY  RATE  OBSTACLE  150,000 

PROPERTY  DAMAGE  RATE  OBSTACLE  15,000 

DEATH  RATE  FOR  OBSTACLES  1000,000 

ACCIUENT  RELATED  COST  DATA 

C03T  PER  DEATH  -50000,000 

COST  OF  PtRSONAL  INJURY  2000,000 

PROPERTY  DAMAGE  <*00.000 
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The  total  economic  response  is  then  computed  and  the  results  of 
the  risk  analysis  for  each  of  the  bridge  designs  is  shown  in  Table  19. 
The  300-foot  bridge  with  a  315-foot  embankment  height  has  a  lower  total 
economic  response  than  the  present  bridge  design,  the  500-foot  bridge  with 
a  315-foot  embankment  height.  This  is  due  to  the  fact  that  the  savings 
in  bridge  cost  of  the  300-foot  bridge  over  the  500-foot  bridge  is  greater 
than  the  increase  in  risk.  This  can  be  readily  seen  by  examining  the 
risk  and  bridge  construction  cost  columns.  For  a  300-foot  bridge  with 
an  embankment  of  315  feet,  the  risk  component  of  the  total  economic 
response  is  approximately  $535  greater  than  the  risk  component  for  the 
present  bridge  configuration.  On  the  other  hand,  a  savings  of  over  $20,000 
per  year  in  construction  cost  is  realized  by  building  the  300-foot  bridge 
rather  than  the  present  bridge  design.  The  data  used  in  this  example 
problem  is  representative  and  not  exact  for  the  given  Tallahalla  Creek 
bridge  site  and  therefore  we  are  not  stating  specifically  that  the  300-foot 
bridge  length  design  is  better  than  the  present  design  although  this  would 
appear  to  be  the  case  based  on  the  representative  data  set. 

The  plot  option  of  the  RISK  module  was  used  to  generate  three 

plots  of  the  information  shown  in  Table  19.  In  the  first  plot  (Figurel8) 

-3        -3 
the  total  economic  response  times  10   (TER  *  10  )  for  the  500-foot 

bridge  category  is  plotted  versus  embankment  height.  The  315-foot 

embankment  height  gives  a  lower  response  than  the  310-foot  and  320-foot 

embankment  heights  for  a  500-foot  bridge.  The  economic  response  of  the 

315-foot  embankment  height  is  lower  than  that  of  the  310-foot  because 

the  savings  in  construction  cost  of  the  310-foot  embankment  height  over 

the  315-foot  configuration  ($5,400)  fails  to  offset  the  larger  risk 

component  of  the  310-foot  embankment  height  configuration  ($27,000). 

The  315-foot  embankment  height  economic  response  is  lower  than  that  of 

the  320-foot  because  increase  in  construction  cost  of  the  320-foot 

configuration  over  the  315-foot  configuration  ($6,000)  is  larger  than 

the  savings  of  the  320-foot  configuration  risk  over  the  risk  of  the 


86 


Table  19.  Results  of  Economic  Analysis 


RJS*    FACTOR    ANALYSIS    RESULTS 

BHlCliE    CONFIGURATION    (*T) 
tE  Mi  t  H  ErttKT,    HT, 


JuO.OuO 

500,000 

500,000 

900,000 

1000,000 

1  0  y  0  ,  0  y  0 

1003,000 


515,000 
319,000 

J10.000 
320,000 
315,000 
310,000 
320, COO 


ECONOMIC  RESPONSE  COMPONENTS  (0OLLAR8) 
CONSTRUCTION  COST      RISK      TyTit  ecONUMlC  HCSPONSE 


92172,305 

7271S.814 

67312.107 

787S0  .9S"* 

122919, 755 

119293,935 

124963,337 


1483,047 
948,476 

33251.00* 
606,000 
658,443 

28491.836 
495,000 


53655,352 
73664,292 

100563.107 

79356,984 
123578,197 
147785,771 
127478,337 


Hl3<  BASED  ON  PyOOO  EVENTS 
flow      VULUit(CFS)      RECU«HEnC£  INTERVAL 

1  lo200,000  ,010 

2  19TJ0.000  ,020 

3  23500,000  ,010 
«        27500,000  ,010 


CONST9JCTIUN  COST  CALCULATED  AS  EUUIVALENT  ANNUAL  COST 
BASED  ON  PHOJECT  H^'fe    iOO.OOOVEAKS      INTEREST  HATE      7,000 


PERCENT 


TOTAL  tCONUilC  RtSPONSt  EQUALS  SOU  OF 
CONSTRUCTION  CUST  ANy  KI3K 


Calculation  based  on  incomplete  information 
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315-foot  configuration  ($340).  In  the  second  plot  (Figure  19)  the  total 

-3         -3 
economic  response  time  10   (TER  *  10  )  for  the  1000-foot  bridge 

category  is  plotted  versus  embankment  height.  The  results  are  similar 

to  those  discussed  for  the  500-foot  bridge  category.  In  the  third 

-3         -3 
plot  (Figure  20)  the  total  economic  response  times  10   (TER  *  10  )  is 

plotted  versus  bridge  length  for  all  bridge  configurations  having  an 

embankment  height  of  315  feet.  The  increase  in  total  economic  response 

with  increasing  bridge  length  at  the  315-foot  embankment  height  is 

linear.  This  is  due  to  the  fact  that  the  construction  cost  components 

of  the  total  economic  response  dominate  the  risk  components  for  all 

configurations  plotted,  and  the  construction  cost  is  linearly  monotonic, 

increasing  with  bridge  length  for  a  fixed  embankment  height. 

In  this  example,  the  dominance  of  the  construction  cost  components 
over  the  risk  components  for  all  bridge  designs  considered  is  due  to 
the  bridge  site  location  in  a  rural  area  where  the  backwater  damages 
are  minor.  In  the  following  chapter  on  Sensitivity  Analysis,  a  site 
will  be  considered  where  the  backwater  damage  is  great  and  the  risk 
component  of  the  economic  response  is  comparable  to  the  bridge  construction 
cost. 
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V.  SENSITIVITY  ANALYSIS 


In  the  development  of  the  methodology  for  performing  risk  analysis 
of  highway  stream  crossings,  assumptions  were  made  for  several  physical 
processes.  For  example,  estimates  were  made  from  available  information 
for  such  design  parameters  as  traffic  level  and  characteristics  >  and 
hydrologic  data.  These  assumptions  and  estimates  are  subject  to  errors. 
Hence  the  use  of  this  information  in  calculating  the  total  economic 
response  at  a  bridge  site  may  induce  solution  errors,  which  in  turn 
may  affect  the  decision-making  process. 

The  number  of  parameters  involved  in  the  risk  analysis  is  large 
and  the  governing  hydrologic  and  hydraulic  processes  are  complex.  More- 
over, the  economic  data  used  for  the  analysis  varies  from  site  to  site 
and  also  with  time.  The  uncertainties  associated  with  the  risk  analysis 
need  not  discourage  the  use  of  the  method  as  long  as  the  limitations 
inherent  in  the  method  of  analysis  are  identified.  Sensitivity  analysis 
is  the  means  by  which  these  limitations  may  be  assessed.  The  results 
of  sensitivity  analysis  permit  the  designer  to  deal  with  the  uncertainty 
presented  in  the  method  of  analysis.  Sensitivity  analysis  identifies 
critical  variables  and  factors  for  the  purpose  of  singling  them  out 
for  careful  consideration  in  the  design  process  to  reduce  uncertainty. 
It  also  serves  the  purpose  of  establishing  research  priorities  and 
resolves  estimating  problems  associated  with  sensitive  design  factors. 

Virtually  every  parameter  and  measurement  associated  with  bridge 
design  is  subject  to  some  degree  of  estimation.  Basically  there  are  three 
major  groups  of  parameters  governing  the  total  cost  for  a  crossing  site. 
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These  are:  (1)  hydrologic  -  hydraulic,  (2)  engineering,  and  (3)  economic 
variables.  In  .a  previous  study  by  WRE  (13)  for  22  culvert  sites,  it 
was  found  that  six  of  the  factors  in  the  three  groups  appear  to  be 
relatively  significant  in  the  economic  response.  These  factors  are: 

•  interest  rate, 

•  unit  cost, 

•  stage-damage  curve, 

•  flow  hydrographs,  and 

•  embankment  erosion. 

This  list  is  the  basis  for  the  selection  of  factors  in  the  present  sensitivity 
analysis.  The  major  parameters  subject  to  analysis  are: 

t  interest  rate, 

•  bridge  cost, 

•  stage-damage  curve  (backwater  damages),  and 

•  hydrology,  i.e.,  flood-frequency  curve. 

Minor  parameters  analyzed  include  (1)  cost  of  time,  and  (2)  average 
daily  traffic  (ADT).  The  results  of  the  sensitivity  computer  runs  are 
shown  in  Table  20.  Run  1  represents  the  baseline  values  (or  the  optimal 
design)  used  in  the  example  problem  of  Chapter  IV.  Runs  2  through  13 
comprise  the  sensitivity  analysis  runs  made  in  this  study.  The  optimal 
solution  remained  the  300-foot  bridge  with  a  315-foot  embankment  height 
for  all  the  sensitivity  analysis  runs.  The  effect  of  changes  in  the 
sensitivity  analysis  parameters  on  the  total  economic  response  of  this 
optimal  bridge  design  is  presented  in  Figure  21  .  The  relative  effect  of 
each  parameter  varied  in  the  sensitivity  analysis  on  the  total  economic 
response  of  the  optimal  bridge  design  is  evaluated  in  Table  21  .  In 
Table  22  the  sensitivity  analysis  parameters  are  ranked  by  order  of  impact 
on  the  total  economic  response  of  the  optimal  bridge  design. 
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Table  21.  Effect  of  Parameter  Changes 
on  Total  Economic  Response 


Run 


Parameter 


ATER/TER  AP/P 


Avg.  Relative 
Effect  on 
Response 


2  Interest  Rate  5% 

3  Interest  Rate  9% 

4  Cost  of  Time  ($3.64/hr) 

5  Cost  of  Time  ($4.57/hr) 

6  Average  Daily  Traffic  Level 
(12000  veh/day) 

7  Average  Daily  Traffic  Level 
(20000  veh/day) 

8  +25%  Flow  for  Given  Frequency 

9  -25%  Flow  for  Given  Frequency 

10  +25%  Bridge  Cost 

11  -25%  Bridge  Cost 

12  Baseline  Backwater  Damages  X  10 

13  Baseline  Backwater  Damages  X  100 


0.0127 


0.1205 

0.072 

1.178 

1.178 

0.0115 

0.0127 


0.0127 


0.096 
1.178 
0.0121 


Relative  Effect  on  Response 


Table  22.  Ranking  of  the  Parameters  of  the  Sensitivity 
Analysis  by  Impact  on  Total  Economic  Response 


Parameter 


Rank 


Bridge  Cost  1 

Interest  Rate  2 

Flow  Frequency  Curve  (Gumbel  Plot)     3 

Average  Daily  Traffic  Level  4 

Flood  Losses  5 

Cost  of  Time  6 


96 


Bridge  cost  has  the  largest  impact  on  the  total  economic  response 
of  the  optimal  crossing  design.  To  produce  the  +25  percent  changes  in 
bridge  cost,  the  bridge  regression  coefficients  B(l)  and  B(6),  the  cost  of 
roadway  per  foot  (RWCST),  and  the  embankment  fill  cost  tCY)  are  all 
increased  and  decreased  by  25  percent.  The  cost  of  roadway  per  foot  and 
the  embankment  fill  cost  are  also  used  in  determining  erosion  losses  and 
thus  have  an  effect  on  the  risk  component  of  the  total  economic  response. 
A  one  percent  change  in  B(l),  B(6),  RWCST  and  CY  produces  a  1.178  percent 
change  in  the  total  economic  response  of  the  optimal  solution. 

The  interest  rate  used  in  determining  the  equivalent  annual  bridge 
construction  cost  is  the  next  most  critical  parameter.  The  equivalent 
annual  cost  is  determined  by  multiplying  the  bridge  cost  by  the  capital 
recovery  factor.  The  total  economic  response  is  the  sum  of  the  equivalent 
annual  cost  of  the  bridge  and  the  risk.  A  one  percent  change  in  interest 
rate  produces  a  0.987  percent  change  in  the  total  economic  response  of 
the  optimal  solution. 

The  nature  of  the  flow  frequency  curve  has  a  significant  impact 
on  the  total  economic  response  of  the  optimal  bridge  design.  The  flow 
frequency  curves  used  in  the  sensitivity  analysis  are  presented  as  the 
Gumbel  plots  shown  in  Figure  22.  For  a  given  recurrence  interval  the 
flow  is  increased  and  decreased  by  25  percent  in  order  to  develop  the  plus 
and  minus  25  percent  straight  lines  shown  in  Figure  22.  The  return  weights 
used  to  develop  the  risk  component  of  the  total  economic  response  are 
determined  from  the  Gumbel  plot  of  the  flood  frequency  information.  The 
evaluation  of  the  return  weights  is  also  based  on  the  flood  events 
considered  in  the  risk  analysis  and  is  shown  in  Table  23.  A  one  percent 
change  in  the  Gumbel  plot  of  Figure  22  produces  a  1.096  percent  change 
in  the  total  economic  response  of  the  optimal  bridge  design. 
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Table  23.  Evaluation  of  Return  Weights  for 

Flood-Frequency  Sensitivity  Analysis 


Flow  (cfs) 


Loss  for 
Flow  Q  ($) 


Return 
Interval 

(yr) 


1=1, ..5    i=l,..4      Return 

p[Q>QJ       p[Qi<Q<Qi+1]     ^ight 


CD 


15000 
16200 
14700 


S   23500 


CO 


27500 


15000 


"16200 
•14700 
"23500 
"27500 


20 

25 

50 

100 

200 


0.05 
0.04 
0.02 
0.01 
0.005 


0.01 
0.02 
0.01 
0.005 


0.01 
0.02 
0.01 
1—0.01 


-PfCHjJ  =  0.005- 


Risk  =  0.01  L1620Q  +  0.02  L147Q0  +  0.01  I  +  0.01  L27500 


Baseline  Return  Weight  Set  = 


0.01 
0.02 
0.01 
0.01 
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Table  23.  Evaluation  of  Return  Weights  for  Flood- 
Frequency  Sensitivity  Analysis  (Cont.) 


Flow  (cfs) 


Loss  for 
Flow  Q  ($) 


Return 
Interval 

(yr) 


i=1, ..5    i=l,..4      Return 
P[Q>QJ    P[Q1<Q<Q1+1]   Weight 


CD 
O 

s- 

0) 
Q_ 

LO 
C\J 

+ 


15000 
16200 
19700 
23500 
27500 


■15000 


16200 


19700 


23500 


L27500 


10 

12.5 

20 

35 

70 


0.1 

0.08 

0.05 

0.028 

0.014 


0.02 
0.03 


>- 

} 0^022 

\ 0.014- 


0.02 
0.03 
0.022 
—  0.028 


P[Q>Q5]  =  0.014-1 


Risk  =  0.02  L1620Q  +  0.03  L1970Q  +  0.022  I  +  0.028  L2750Q 


+25%  Return  Weight  Set 
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Table  23.     Evaluation  of  Return  Weights  for  Flood- 
Frequency  Sensitivity  Analysis   (Cont.) 


Q 
Flow  (cfs) 

Loss  for 
Flow  Q  ($) 

Return 
Interval 

(yr) 

15000 

L15000 

50 

!  16200 

i 

L16200 

70 

;  19700 

L19700 

190 

!  23500 

L23500 

500 

27500 

Lo-7c:nn 

1700 

i=l,.._5  i=l,.. 4  Return 

p  B>QiJ         p[Qi<Q<Qi+1]       ^ight 


CD 
O 

S- 

<u 

Q- 

I 


0.02 

0.014 

0.0053 

0.002 

0.0006 


0.006 
■0.0087 
-0.0033 
■0.0014- 


LP[Q>Q5"]   =  0.0006- 


0.006 
0.0087 
0.0033 
0.002 


Risk  =  0.006  L1620Q  +  0.03  L19700  +  0.022  L235()0  +  0.028  L27500 


-25%  Return  Weight  Set  = 


0 

006 

0 

0087 

0 

0033 

0 

002 
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The  average  daily  traffic  level  in  vehicles  per  day  is  used  to 
determine  the  total  dollar  amount  of  traffic-related  losses  caused  by 
inundation  of  candidate  bridge  designs.  The  total  dollar  amount  of  traffic- 
related  losses  forms  one  component  of  the  total  risk  which  directly  affects 
the  total  economic  response.  A  one  percent  change  in  average  daily  traffic 
level  is  small  at  0.0127  percent. 

The  losses  due  to  backwater  created  by  flow  constriction  of  the 
bridge  are  a  component  of  the  total  risk.  The  baseline  backwater  damages 
are  in  the  order  of  103  dollars  in  magnitude  and  are  those  associated 
with  a  real  site  such  as  Tallahalla  Creek.  In  the  sensitivity  analysis 
these  backwater  damages  are  first  multiplied  by  a  factor  of  10  to 
represent  an  urban  flood  plain  damage  pattern  and  then  by  a  factor  of  100 
to  represent  industrial  flood  plain  damage.  The  total  economic  response  of 
the  optimal  bridge  design  is  quite  insensitive  to  changes  in  backwater 
damages  even  in  the  case  of  industrial  flood  plain  development  for  which 
the  risk  components  of  the  total  economic  response  are  on  the  same  order 
of  magnitude  as  the  bridge  construction  cost.  For  a  one  percent  change 
in  backwater  damages  a  change  of  0.0121  percent  in  the  total  economic 
response  of  the  optimal  solution  is  observed. 

The  cost  of  time  in  dollars  per  hour,  like  the  average  daily 
traffic  level,  is  used  in  determining  the  total  dollar  amount  of  traffic 
delay  losses  incurred  by  overtopping  of  the  bridge  embankments  and 
subsequent  embankment  erosion.  The  total  dollar  amount  of  traffic  delay 
losses  comprises  one  component  of  the  total  risk  wfrich  in  turn  ts  used 
to  evaluate  the  total  economic  response.  The  input  of  a  one  percent 
change  in  the  cost  of  time  in  dollars  per  hour  on  the  total  economic 
response  of  the  optimal  bridge  design  is  quite  small  at  0.0013  percent. 

In  conclusion,  the  optimal  bridge  design  for  the  Tallahalla  site 
is  the  300-foot  bridge,  315-foot  embankment  height  configuration  for  all 
changes  in  the  sensitivity  analysis  parameters.  Interest  rate  and  bridge 
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cost  have  a  major  influence  on  the  total  economic  response  of  the  optimal 
bridge  design.  Changes  in  the  flood  frequency  curve  significantly  impact 
the  total  economic  response  of  the  optimal  solution,  while  changes  in 
average  daily  traffic  level,  backwater  damages,  and  cost  of  time  produce 
only  minor  changes  in  the  total  economic  response. 
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VI.  SUMMARY,  CONCLUSIONS  AND  RECOMMENDATIONS 


SUMMARY 

An  engineering  systems  analysis  technique  has  been  developed  to 
facilitate  the  economic  risk  analysis  for  highway  stream  crossings.  The 
technique  is  based  on  a  probabilistic  approach,  incorporating  hydrologic, 
hydraulic  and  economic  factors  in  the  engineering  analysis  of  bridge 
waterways.  It  relates  the  initial  construction  costs,  flood-related 
highway  losses  and  damage  to  properties  from  backwater  to  crossing  design 
and  flood  frequency. 

The  basic  concept  of  the  risk  analysis  of  highway  stream  crossing 
is  a  method  to  determine  the  economic  measurement  of  various  crossing 
designs  at  a  bridge  site.  The  economic  measurement  used  in  this  study 
is  the  total  economic  response,  defined  as  the  sum  of  project  construction 
cost  and  risk,  on  an  annual  basis.  A  design  scheme  giving  the  least 
total  economic  response  is  the  optimal  design. 

The  annual  project  construction  cost  is  derived  from  the  cost  of 
bridge,  approach  embankments  and  the  pavement,  amortized  over  the 
project. life  at  a  given  interest  rate.  Such  a  cost  is  generally  governed 
by  the  type  of  bridge,  bridge  length,  bridge  width,  clearance,  span, 
geology  or  underlying  river  bed,  height  of  the  approach  embankments,  type 
of  highway,  and  width  of  the  valley  at  the  crossing  site.  Risk  is  the 
sum  of  the  products  of  the  weighted  probability  of  flood  occurrence  and 
flood-related  economic  losses.  The  determination  of  risk  requires  hydrologic 
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hydraulic  and  economic  data  at  and  near  the  highway  crossing  site. 
Tangible  losses  to  be  considered  include  structure  damages,  traffic-related 
(e.g.,  traffic  stoppage  or  delay)  losses,  and  flood  losses. 

Hydrologic  data  are  obtained  from  the  U.  S.  Geological  Survey 
stream  gaging  records.  Flood  data  are  plotted  on  Gumbel  probability  paper 
to  define  the  flood-frequency  curve  at  the  crossing  site.  Hydraulic  data 
such  as  the  water  surface  elevations  and  the  velocities  at  and  near  the 
bridge  site  are  computed  by  a  two-dimensional  Finite  Element  Model. 

The  decision  variables  considered  for  various  design  alternatives 
in  this  study  are  limited  to  the  bridge  length  and  embankment  height. 
Other  minor  variables  such  as  bridge  type  and  crossing  type  may  be  incor- 
porated in  the  analysis  when  situations  warrant.  The  complexity  of  the 
decision  matrix  increases  as  the  number  of  decision  variables  increases. 

The  method  of  risk  analysis  is  implemented  by  a  model  in  the  form 
of  a  computer  program  using  FORTRAN  program  language.  Details  of  the 
model  description  and  its  input  specifications  are  presented  in 
Appendix  B. 


CONCLUSIONS 

The  following  conclusions  have  been  drawn  from  the  present  study: 

1.  The  study  has  illustrated  that  the  application  of  systems 
analysis  techniques  to  provide  the  means  for  design 
engineer's  decision-making  process  in  highway  stream  crossings 
is  feasible  and  practical.  This  report  has  successfully 
demonstrated  a  methodology  to  deal  with  a  complex  problem 
by  a  systematic  simplification  of  the  analysis. 
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2.  Results  of  sensitivity  analysis  for  the  parameters 
tested  reveal  that  the  relative  importance  in  design 
parameters  affecting  the  total  economic  response  for 
the  example  site  is  in  the  following:  bridge  cost, 
interest  rate,  flow  frequency  curve,  ADT,  flood  losses 
and  the  cost  of  time.  For  bridges  subject  to  inundation, 
results  of  sensitivity  analysis  may  vary. 

3.  Because  the  example  site  is  located  in  a  rural  area,  the 
results  show  that  the  bridge  cost  outweighs  the  risk. 
Such  an  observation  may  not  be  true  for  bridge  sites 
located  in  urban  areas. 

4.  With  the  limited  data  available  at  the  Tall  anal  la 
Creek  site,  the  results  of  the  example  problem  have 
indicated  that  the  present  method  gives  a  more 
economic  design  than  the  conventional  method.  This 
is  consistent  with  the  results  arrived  at  in  the 
previous  study  for  box  culverts  (13). 


RECOMMENDATIONS 

The  research  effort  described  herein  encountered  numerous  problems, 
the  solution  to  which  should  be  pursued  simultaneously  on  a  twofold  front.  The 
estimation  of  such  important  variables  as  pier,  abutment,  and  channel  scour, 
embankment  erosion,  and  bridge  failure  mechanisms,  require  work  on  both  simplistic 
empirical  formulations  and  basic  research  to  develop  methodology  that  is 
generally  applicable  to  the  wide  range  of  conditions  normally  encountered.  Of 
course,  this  should  not  delay  the  application  of  current  technology  to  solving 
present  day  problems.  It  is  only  in  this  way  that  the  needed  data-base  can  be 
assembled  from  which  new  techniques  can  be  established  and  tested.  In  most 
cases  the  methodology  employed  in  this  report  is  not  generally  applicable  to  all 
cases  but  must  be  calibrated  by  a  trial  and  error  process  to  the  bridge  site  under 
study.  Ideally,  as  more  information  is  gained,  the  difficulty  in  estimating 
important  parameters  will  be  alleviated  by  a  compilation  of  a  set  of  "Handbook" 
values  applicable  to  a  wide  range  of  situations.  Specific  recommendations  in 
this  regard  follow: 
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1.  A  more  phenomenologically  oriented  method  for  calculating 
embankment  erosion  is  needed.  The  method  used  here  is  a  log- 
linear  relationship  developed  as  a  function  of  velocity. 
Additional  study  is  required  to  either  improve  this  empirical 
relationship  or  to  develop  alternative  approaches. 

2.  Presently,  only  pier  scour  is  included  in  the  Risk  analysis. 
Abutment  and  channel  scour  should  also  be  included  in  the  analysis. 
An  approach  such  as  proposed  by  Laursen  (15)  would  be  a  useful 
addition  to  the  analysis. 

3.  Bridge  damage,  caused  by  inundation,  is  handled  in  only  a 
superficial  manner.  However,  the  method  accurately  quantifies 
the  amount  of  information  presently  available  describing  this 
problem.  This  presents  an  area  where  basic  research  and  evaluation 
of  data  is  urgently  needed. 

4.  In  order  to  move  the  model  from  a  research  tool  to  a  field 
implementable  one,  it  is  recommended  that  additional  sites  be 
analyzed  with  the  present  techniques.  This  will  provide  a  basis 
for  the  improvement  of  present  methods  and  will  allow  valuable 
input  by  field  engineers  into  the  structure  of  the  model. 

5.  The  hydraulic  model  should  be  simplified  in  order  to  reduce  the 
cost  of  analysis.  Such  a  simplification  could  be  made  by 
compiling  hydraulic  data  covering  a  wide  range  of  field  design 
conditions  by  the  use  of  the  Finite  Element  Model  developed  under 
this  contract. 
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GLOSSARY  OF  MATHEMATICAL  TERMS 


a„  =  Constant 
o 

a-.  =  Coefficient  for  area  (x-, ) 

a~  =  Coefficient  for  area  (x2) 

a3  =  Coefficient  for  area  (x3) 

a«  =  Coefficient  for  area  (x*) 

a5  =  Coefficient  for  area  (x5) 

A  =  Intercept  coefficient 
o        r 

A-,  =  Slope  coefficient 

b  =  Eroded  roadway  width 

b  =  Width  of  retangular  pier 

B  =  Embankment  width 

B3  =  Abutment,  footing  and  pier  costs 

B*  =  Protective  measure  costs  (spur  dikes  and  riprap) 

C  =  Cost  adjustment  factor 
a 

C.p  =  Original  cost  of  fill 

Cr  =  Original  cost  of  roadway 

C  =  Cost  of  bridge  superstructure 

CRF  =  Capital  recovery  factor 

D.  =  Loss  incurred  due  to  occurrence  of  state  of  nature  j 

•J 

d  =  Flood  head  above  roadway 

df  =  Pile  length  required  for  nominal  scour  depth 
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d  =  Depth  of  scour 

E  =  Embankment  erosion  (tons/ft/day) 

G  =  Cost  per  unit  pile  length 

H  =  Fill  height 

H.  =  Fill  height  at  station  i 

i  =  Interest  rate 

L  =  Bridge  length 

L.  =  Length  of  bridge  embankment  element  i 

L.  =  Bridge  damage 

P.  =  Probability  of  occurrence  of  state  of  nature  j 

Q  =  Flow 

Q  =  Peak  flow 

r.  =  Recurrence  interval  associated  with  given  flow  Q 
l  3        ^r 

R.  =  Risk  for  strategy  i 
RW  =  Roadway  length 

R  =  Risk 
S-, ,  S«  =  Side  slope  of  embankment  fill 

T  =  Time  of  overtopping 

T.  =  Flood  duration 

T  =  Time  of  concentration 
c 

T  =  Time  to  peak 

T  =  Time  between  peak  flow  and  end  of  hydrograph 
TER  =  Total  economic  response  function 
V  =  Velocity  in  ft/sec 
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V  =  Erosion  threshold  velocity 
y  =  Submergence  of  bridge  deck 


a  =  Empirical  constant  for  embankment  erosion  for  cohesive  soil 
3  =  Empirical  constant  for  embankment  erosion  for  cohesionless  soil 
Y  =  Correlation  coefficient 

t,   =  Coefficient  for  bridge  damage  due  to  inundation 
a      =   Standard  error  for  bridge  cost  regression 
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APPENDIX  A  FEDERAL  HIGHWAY  ADMINISTRATION  POLICY  FOR  HYDRAULIC  DESIGN  OF 
HIGHWAY  ENCROACHMENTS  ON  FLOOD  PLAINS 
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VOLUME 
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CHAPTER 

7 

BRIDGE,    STRUCTURES   AND    HYDRAULICS 

SECTION 

3 

HYDRAULICS,    EROSION   CONTROL   AND   WATER   QUALITY 

Subsection  2 


HYDRAULIC  DESIGN  OF  HIGHWAY  ENCROACHMENTS  ON 
FLOOD  PLAINS 


Transmittal    30 

Par.    1. 

Purpose                                                               May    29,    1974 

2. 

Authority                                                          HNG-31 

3. 

Policy 

4. 

Definitions 

5. 

Design   Standards 

6. 

Studies    and   Reports 

7. 

Federal    Participation    in   Construction   Costs 

1.      PURPOSE 

The   purpose    of   this    directive    is    to   prescribe   policies    and 
procedures    for   hydraulic   designs    for   highway    projects 
constructed  with   Federal-aid   funds    and  projects    under    the 
direct    supervision   of   the    Federal    Highway   Administration, 

AUTHORITY 

This  is  issued  under  authority  of  23  U.S.C.  109(a),  315, 
23  CFR  1.32,  49  CFR  1.48(b)(8)  and  implements  Executive 
Order  11296  (August  10,  1966,  3  CFR,  Part  II). 


3.   POLICY 


Pursuant    to    Executive    Order    11296,    it    is    the   policy    of 
the    Federal    Highway   Administration    to    encourage    a   broad 
and   unified   effort    to   prevent    uneconomic ,    hazardous    or 
unnecessary    use    and   development    of   the    Nation' s    flood 
plains ,    and   in   particular    to    lessen    the    risk   of  flood 
losses    in    connection   with   Federally- financed   or    supported 
improvements ;    and    to    comply   with    the    "Flood   Hazard 
Evaluation    Guidelines   for   Federal    Executive   Agencies," 
May    1972,    published   by    the    Water   Resources    Council, 
(See  Attachment  1) 
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b.  It    is    the    policy    of   the    Federal    Highway   Administration 
that,    where    practicable ,    highway    locations    shall 
avoid    areas    subject    to    flooding . 

4.  DEFINITIONS 

a.  The    term    "basic    flood"    shall    mean    the    1 00-year 
floodo 

b.  The    term    "conveyance    of   the    basic    flood"    shall    mean    the 
ability    to    accommodate    passage    of   the    1 0 0-year   flood, 

(1)  Conveyance    may    be    through    structures    or   both 
through    structures    and   over    the    highway . 

(2)  Conveyance    along    a    highway    may    include    inundation 
of    the    highway . 

c.  The    term    "design    flood"   shall   mean    the    peak    discharge, 
volume    (if  appropriate ) 3    and   stage    or   wave    crest 
elevation    of   the    flood   associated   with    the    recurrency 
interval    selected   for    the    design    of  a    highway 
encroachment    on    a    flood   plain.       By    definition,    the 
highway    will    not    be    subjected    to    inundation   from    the 
stage  of   the    design   flood. 

d.  The    term    "flood   plain"    shall    mean    (1)    the    valley    area 
adjacent    to    a    stream   or   river   subject    to    overflow , 
02)    an   area    adjacent    to    a    lake,    an    estuary ,    an    ocean, 
or    similar   body    of  water    subject    to    high    tides,    surges, 
tsunamis ,    or   any    combination    of   these;    or    (Z)    an 

area   where    the    path    of   the    next   flood   flow    is 
unpredictable ,    as    in    a   debris    cone,    an    alluvial    cone 
or   fan,    a    debris    slope,    or   a    talus. 

e.  The    term    "encroachment"   shall    mean    a    highway    and/or 
appurtenant   feature    within    the    limits    of  a   flood  plain. 

5.  DESIGN  STANDARDS 

a.  Where    encroachment    on   a   flood  plain    is    necessary ,    an 
evaluation,    using    the    basic   flood,    shall    be   made    of 
the    flood   hazard    to    the    highway ,    and   the    effect    of 
the   proposed   highway    on    the   flood   hazard   to    other 
property ,    stream   stability ,    and    the    stream   and  flood 
plain    environment .       (See  Attachment  2) 

b.  The  design  flood  selected  for  each  encroachment  and  the 
design  for  conveyance  of  the  basic  flood  should  be 
supported  by  the  following,  where  applicable:   an 
incremental  analysis  of  estimated  construction  costs; 
probable  property  damage  including  damage  to  the 
highway;  the  cost  of  traffic  delays;  the  availability 
of  alternate  routes,  emergency  supply  and  evacuation 
routes;  and  with  consideration  of  the  potential  for 
loss  of  life  and  of  budgetary  constraints. 
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k.  Buildings 3    waste    treatment    facilities    and    other    high, 
cost    installations    or   potential    major    source    of 
pollution    for    highway    rest    areas    shall    be    located 
outside    the    flood   plain    or   flood   proofed   against 
damage    from    the    basic   flood. 

1.  Where    highway    fills    are    to    be    used    to    form    dams    or 
levees3    appropriate    permits  3    with    regard    to    water 
rights,    impoundments 3    and   diversions    shall    be    obtained. 
Prior    to    authorization    by    the    Division    Engineer 3     the 
hydrologic  3    hydraulic  3    and    structural    design    of   the 
fill    and   appurtenant    spillways ,  shall    have    the    approval 
of    the    State    or   Federal    agency    responsible    for    the 
safety    of  dams    or    like    structures    within    the    State. 

6.  STUDIES  AND  REPORTS 

a.  Reports    shall    be   prepared   by    the    State    highway 

department    or    its    agent    showing    (1)    hydrologic    and 
hydraulic    data    and   design    computations 3     (2)    the    analysis 
of   the    highway    effect    on    stream    stability  3    and    (3) 
the    favorable    or    adverse    effects    on    the    stream 
environment. 

be  The    reports    and   design    computations    shall    be    retained 
in    the    highway    agency's    permanent    design    files    for 
the    project3    and   upon    request    made    available    to    the 
Federal    Highway   Administration   for   review.       Copies 
of   these    reports    shall    accompany    preliminary    plans 
which    require    regional    or    Washington    office    approval . 

c.  The   reports    shall    be    commensurate   with    the    importance 
of  the    structure    and    the    risk    to    life    and  property . 
Less    comprehensive    reports    are    appropriate    for 
encroachments    on   minor   streams    and  at    locations    where 
the   risk   of  property    damage    or    loss    of   the    highway 
is    small. 

7.  FEDERAL  PARTICIPATION  IN  CONSTRUCTION  COSTS 

a.  Federal-aid  participation    in   cooperative   projects 
using    the    highway   fill    to   provide   for   flood   control 
or   water   resource    development   will    be    limited   to 
the   pro   rata   share    of   the    estimated   cost    to    construct 
the    highway    only. 
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c.  All    highways    that    encroach    on    flood   plains,    bodies    of 
water    or    streams    shall    be    designed    to    permit    conveyance 
of   the    basic    flood   without    causing    significant    damage 
to    the    highway ,     the    stream,    body    of   water,    or    other 

property , 

d.  *  At    each    location    where    the    highway    will    encroach    upon 

a    flood   plain,    project    plans    shall    show    (1)    the 
magnitude ,    frequency ,    and   water    surface    elevations    for 
the    design    flood    and    the    basic    flood,     if   different 
from    the    design    flood,    and    (2)    the    magnitude ,    water 
surface    elevations    and    date    of   occurrence    of   the 
flood    of   record,     if   greater    than    the    basic    flood. 

e.  The    effects    of   existing    flood    control    channels ,     levees, 
and    reservoirs    shall    be    considered    in    estimating    the 
peak   discharge    and   stage    for   all    floods    considered   in 
the    design. 

£.  All    highway s    on    the    Interstate    System    that    encroach   on 
flood   plains    shall    be    designed    to    avoid    inundation    of 
the    highway    from   floods    at    least    as    great    as    the    50- 
year    flood. 

g.   Encroachments  on  flood  plains  by  ramps  for  traffic 
interchange  with  Interstate  highways  may  be  designed 
to''  permit  inundation  of  the  ramp  from  floods  smaller  than 
the  50-year  flood  where  traffic  service  would  not  be 
unduly  reduced  and  where  the  highest  standards  of 
service  are  not  warranted  for  the  crossroad. 

h.   Encroachments  on  flood  plains  by  Interstate  highway 
service  roads  or  frontage  roads  should  be  designed  to 
permit  inundation  of  the  road  from  floods  smaller  than 
the  50-year  flood  except  at  locations  where  lower 
design  standards  for  the  service  road  would  interfere 
with  the  performance  of  hydraulic  structures  for  the 
Interstate  roadway  or  where  higher  standards  can  be 
shown  to  be  warranted. 

i.   Highways,  other  than  those  on  the  Interstate  System, 
which  encroach  on  flood  plains,  may  be  designed  to 
permit  inundation  by  floods  smaller  than  the  50-year 
flood,  where  lower  standards  are  warranted. 

j.   Ramps,  parking  and  picnic  areas  for  rest  areas  may  be 
designed  to  be  inundated  by  overbank  flow  in  the  flood 
plain,  dependent  upon  the  expected  duration  of  flooding, 
expected  cleanup  and  repair  costs  and  the  hazard  to 
users.   Precautions  should  be  taken  to  insure  that 
benches,  tables,  trash  barrels  and  other  furnishings 
do  not  become  waterborne  debris  during  floods. 
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The  following  extracts  from  "Flood  Hazard  Evaluation 
Guidelines  for  Federal  Executive  Agencies,"  United  States  Water 
Resources  Council,  May,  1972,  are  applicable  to  the  location, 
design  and  construction  of  Federal  and  Federal-aid  highways. 
The  WRC  publication  has  been  distributed  by  various  sources. 
Additional  copies  are  available  from  FHWA  on  request. 

"(1)   Determine  first,  ...,  whether  there  is  any  need  to 
evaluate  the  flood  hazard  at  the  site  or  structure  location  being 
considered. " 

"(3)   Use  the  following  to  identify  and  evaluate  the  flood 
hazard :" 

"(3A)   The  100-year  flood  as  the  basic  flood:" 

"(3B)   The  flood  hazard  zone,  defined  as  the  area  inundated 
by  the  basic  flood:" 

"(3C  2)  In  the  Case  Approach,  a  procedure  to  assure  that  any 
encroachment  on  the  flood  plain  will  permit  conveyance  of  the 
basic  flood  without  increasing  flood  heights  or  velocities  to 
an  extent  which  would  cause  significant  upstream  or  downstream 
damage  to  existing  or  reasonably  anticipated  future  development:" 

"(3D)   Floods  greater  or  less  than  the  basic  flood  as 

appropriate . " 

"(4)   Determine  whether  there  are  existing  laws  or  statutes 
of  the  Federal  Government,  rules  or  regulations  of  other  Federal 
agencies,  or  laws,  statutes,  ordinances,  etc.,  of  State  or  local 
governments  that  provide  standards  for  regulation  of  the  floodplain 
under  study.   In  cases  where  those  standards  are  either  more 
stringent  than  those  based  on  these  guidelines,  or  are  applicable 
to  situations  or  conditions  not  covered  by  these  guidelines,  they 
should  be  considered  for  the  evaluation  of  flood  hazard  in  that 
area ." 

"(5)   Decide  on  the  conditions  under  which  an  evaluation 
must  be  made  to  determine  the  impacts  of  including  or  excluding 
the  use  of  site  in  a  floodplain.   Such  evaluation  must  demonstrate 
clearly  that  the  use  of  the  site  is  to  the  advantage  of  society 
as  well  as  to  the  user  of  such  site." 
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"(6)   Select  the  floods  to  be  used  in  a  flood  hazard  evaluation 
to  fit  conditions  of  the  area  being  investigated." 

"(10)   Determine  the  effects  of  proposed  highway  construction 
in  the  floodplain  and  its  vicinity,  and  of  proposed  upstream  or 
local  flood  prevention  or  control  measures,  if  any,  on  the 
elevations  of  the  evaluation  floods." 

"(12)   Adopt  the  policy  of  discouraging  the  construction  of 
those  roads,  utilities,  and  other  public  facilities  (except  those 
crossing  streams)  within  the  most  hazardous  portions  of  the  flood- 
plain  that  aggravate  flooding  and  encourage  undesirable  developments 
in  that  zone." 

"(14)   Delineate,  or  ensure  the  delineation  of,  on  Federally 

owned  properties,  the  elevation  of  the  100-year  flood,  and  the 

elevations  and  dates  of  occurrence  of  floods  of  record  whose 
magnitudes  should  be  known  by  the  public." 

"(15)   Encourage  State  and  local  agencies  to  keep  a  permanent 
record  of  information  on  each  floodplain  evaluated,  the  flood  hazard 
evaluation  procedure  and  decisions,  the  flood  prevention  or  control 
measures  proposed  for  upstream  or  local  construction,  the  flood 
elevation^delineations . . ." 
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Checklist   for  Drainage   Studies   and  Reports 

(This    is   a  checklist   of  items   to  be   considered   for  inclusion 
in  hydraulic   studies   and   reports.      Numbers   in  parenthesis    indicate 
references    cited   at   the   end  of  the   checklist). 

1.  Preliminary   drainage   surveys (1) 

a.  Investigate  potential   problems 

(1)  Flood  hazard  -    land  use   and  development   in   flood 
plain,    flood  hazard  studies(2) 

(2)  Channel    stability  -   bank   stability,   bends   and  meanders, 
aggradation,    degradation,   necessity   for  channel    change. 

(3)  Effects   on   the   environment,    fish   and  other  wildlife, 
domestic  water  supplies,    recreational    resources,    etc. 

(4)  Debris  and  ice 

(5)  Skew  of  crossing 

b.  Coordination  with  other  agencies 

(1)  Permits  required 

(2)  Existing  and  proposed  water  resource  projects 
(PPM  50-4.2  (Volume  6,  Chapter  1,  Section  1, 
Subsection  4,  of  the  Federal-Aid  Highway  Program 
Manual)) 

(3)  Possibility  of  cooperative  projects 

2.  Design 

a.   Site  data 

(1)   Vicinity  map 

(a)  Purpose  -  to  show  proposed  highway  alinement  and 
reach  of  river,  bends  and  stream  meanders,  general 
flow  directions 

(b)  Type 

1_.   USGS  quadrangle  sheet  or  map  of  equal  detail 
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2_.   Aerial  photo s 

(2)   Site  map: 

(a)  Purpose  -  for  use  in  estimating  flood  flow 
distribution;  to  locate  cross  sections  of  stream; 
to  show  location  of  proposed  encroachment  and 
structure (s) ,  alinement  of  piers,  skew  of  crossing, 
stream  controls,  existing  encroachments,  existing 
highway  structures 

(b)  Type 

1_.   Specially  prepared  map  showing  1-  or  2-foot 
contours,  vegetation  and  manmade  improvements 

2_.      In  some  cases,  cross  sections  normal  to  flood- 
flow  are  acceptable  in  lieu  of  map.   A  minimum 
of  three  sections  are  desirable  including 
cross  sections  upstream,  at  crossing,  and 
downstream 

/  (3)   Existing  structures  (including  relief  or  overflow 
structures) : 

(a)  Locate  existing  structures  with  respect  to  proposed 
crossing  or  encroachment 

(b)  Describe  each  fully,  giving- 

1_.   Type,  including  span  lengths  and  number  of 
spans,  bent  design,  pier  orientation,  culvert 
size,  number  of  cells 

2_.      Foundation  type  (spread  footing,  piling) 
founding  depth 

3_.   Scour  history  at  abutments  and  bents,  culvert 
outlets;  head  cutting;  stream  aggradation, 
degradation 

4_.   Cross  section  beneath  structures,  noting 
clearance  to  superstructure  and  skew  with 
direction  of  current  during  extreme  floods 
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£.   Flood  history,  highwater  marks  (dates  and 
elevation) ,  nature  of  flooding  (including 
overtopping) ,  damages  and  sources  of 
information 

6k   Damage  from  abrasion,  corrosion,  wingwall 
failure,  culvert  end  failure 

(4)   Locate  and  determine  elevations  of  highwater  marks 
along  stream,  giving  dates  of  occurrence 

(a)   Describe  or  list  critical  flood  elevations  of 
interest  in  evaluating  possible  damage  (record 
datum  used) 


(b)   Refer  to  flood  hazard  studies  for  area 


(2) 


(5)  Comment  on  drift,  ice,  nature  of  streambed,  bank 
stability,  bends,  meanders,  vegetative  cover,  land 
use 

(6)  Photographs  showing  existing  structures,  past  floods, 
main  channel  and  flood  plain  to  help  in  evaluating 

a  location  and  documenting  conditions  existing  prior 
to  construction 

(7)  List  factors  affecting  water  stages 

(a)  Highwater  from  other  streams 

(b)  Reservoirs  -  existing  or  proposed  and  approximate 
date  of  construction 

(c)  Flood  control  projects  (give  status) 

(d)  Tides 

(e)  Other  controls 

Hydrologic  analysis  (site  inspection  should  be  made  by 
engineer  making  hydrologic  and  hydraulic  analysis)  (3) (7) 

(1)  List  available  flood  records 

(2)  Determine  drainage  area  above  proposed  construction^-3-' 
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(3)  Evaluate  potential  for  changes  in  watershed  "Characteristics 
which  would  change  magnitude  of  flood  peaks;  e.g., 
urbanization,  channelization 

(4)  Plot  flood- frequency  curve 

(5)  Plot  stage-discharge-frequency  curve 

(6)  Determine  distribution  of  flow  and  velocities  for 
several  discharges  or  stages  in  natural  channel  for 
existing  conditions 

c.   Hydraulic  analysis 

(1)   Design  of  Dridge  waterways  (5j 

(a)  Determine  permissible  upstream  water  surface 
elevations 

(b)  Compute  backwater  for  various  trial  bridge  lengths 
and  discharges 

(c)  Select  design  flood  and  waterway  design 
(Paragraph  4.c.) 

(d)  Provide  for  conveyance  of  100-year  flood 

(e)  Estimate  scour  depth  at  piers  and  abutments (6) 

(f)  Design  riprap  for  bank  protection  and  scour 
attenuation  devices,  if  required. 

(g)  Investigate  need  for  spur  dikes 

(h)   Show  final  layout  in  plan  and  profile 

1_.   Show  design  discharge,  elevations,  and  frequency 
2_.      Show  discharge  and  elevations  of  100-year  flood 

(i)  Comment  on: 

1_.   Selection  of  design  flood 

2.  Conveyance  of  100-year  flood 

3.  Channel  change  (if  provided) 
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4_.   Effects  of  construction 

5_.   Need  for  stream  controls  to  protect  highway 

(2)  Design  of  culverts 

(a)  Determine  allowable  headwater  elevation 

(b)  Compute  and  plot  performance  curves  for  trial 
culvert  sizes(4) 

(c)  Evaluate  abrasion  and  corrosion  potential 

(d)  Select  design  flood  and  culvert  design 
(Paragraph  4.c.) 

(e)  Provide  for  conveyance  of  100-year  flood 

(f)  Evaluate  need  for  debris  control 

(g)  Evaluate  need  for  outlet  protection 

(h)   Investigate  need  for  protection  against  buoyancy 
and/or  failure  by  separation  at  joints (9) 

(i)   Show  final  layout  in  plan  and  profile 

1.  Show  design  discharge,  elevations,  and  frequency 

2.  Show  discharge  and  elevations  of  100-year  flood 
(j)   Comment  on: 

1.   Selection  of  design  flood 
2_.      Conveyance  of  100-year  flood 
3_.   Channel  change  (if  provided) 
4_.   Effects  on  stream  stability 
5_.   Provision  for  fish  passage 

(3)  Design  of  longitudinal  encroachments 
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(a)  Determine   the  effect   of  proposed  encroachment  on 
water-surface  profile  using  various  roadway 
design  alternatives 

(b)  Select   design   flood  and  roadway  design 
(Paragraph  4.c.) 

(c)  Provide  for  conveyance  of  100-year  flood 

(d)  Evaluate   effects   on  scour  and  deposition   in   channel 

(e)  Design  embankment,   bank  and   channel   protection 
needed 

(f)  Show   final    layout   in  plan  and  profile 

1_.      Show  design  discharge,    elevations,    and  frequency 
2_.      Show   discharge   and  elevations  of  100-year   flood 

(g)  Comment  on: 

1_.  Selection  of  design  flood 

2_.  Conveyance  of  100-year  flood 

_3.  Channel  change  (if  provided) 

4_.  Effects  on  stream  stability 

5_.  Effects  on  stream  biology 
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APPENDIX  B 
RISK  MODULE  USER'S  MANUAL 


INTRODUCTION 

A  computer  program  called  RISK  has  been  developed  using  FORTRAN 
programming  language  to  facilitate  the  computation  for  the  risk  analysis 
described  in  this  volume.  The  information  flow  to  the  RISK  module  and 
the  logic  flow  of  the  module  with  subroutine  flowcharts  and  descriptions 
are  discussed  in  the  General  Description  Section.  A  detailed  description 
of  the  necessary  user-supplied  input  data  along  with  a  sample  input  data 
set  are  included  with  a  definition  of  variables  and  a  program  listing. 


GENERAL  DESCRIPTION 

The  RISK  module  accesses  two  types  of  information.  The  first 
is  information  supplied  by  the  FEM  hydraulic  solution  by  means  of  the 
hydraulic  tape  and  the  second  is  the  necessary  user-supplied  data  for 
the  RISK  analysis.  The  information  flow  of  the  RISK  module  is  shown 
in  Figure  23. 

The  RISK  module  consists  of  one  main  calling  program,  four  func- 
tion routines,  and  eight  subroutines  as  indicated  in  Table  24.   The 
detailed  logic  of  the  RISK  program  is  shown  in  Figure  24. 
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MODULE 


USER 
INPUT 


Figure    23       RISK  Module  Information  Flow 
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LOOP  TO   (ao 
FOR  EACH 
DESIGN  SCHEME 


LOOP  TO   (£? 
FOR  EACH 
FLOW  RATE 


START  PROGRAM 
RISK 


CALL 


INRISK 


CALL 


CONCST 


CALL 


RDHYDR 


CALL 


OVTME 


CALL 


EROS  IN 


CALL 


DETOUR 


CALL 


SCOUR 


CALL 


BACKWR 


Q 


0 


SUM  LOSSES 


WEIGHT  LOSSES 


CALL  OUTPUT 


SUM  WEIGHTED 
LOSSES    RISK 


SUM  CONSTRUCTION  COST 
AND  RISK    TOTAL 
ECONOMIC  RESPONSE 


_*_ 


CALL  OUTPT1 


CALL  CONSDA 


CALL  PLOT 


STOP 


Figure  24  Main  Program  RISK 


126 


Table  24   Program  RISK 


Main 


Functions 


Subroutines 


RISK 


DETOUR 
CONCST 
BACKWR 
SCOUR 


INRISK 

RHYDR 

OVTME 

EROSIN 

OUTPUT 

OUTPT1 

CONSDA 

PLOT 


Upon  entrance  to  the  RISK  program  the  INRISK  subroutine  is  called 
to  read  the  user-supplied  input  data.  The  program  flow  enters  the  a2 
loop  and  the  CONCST  routine  is  called  to  determine  the  cost  of  the  given 
bridge  configuration  being  considered.  The  program  then  enters  the 
aj  loop  in  which  a  set  of  flows  is  considered  to  evaluate  the  associated 
risk  components.  For  each  flow  considered,  the  RDHYDR  routine  is  called 
to  read  the  hydraulic  solution  provided  by  the  FEM  module.  Next,  the 
OVTME  routine  is  called  to  calculate  the  overtopping  time  for  each  embank- 
ment section  following  which  the  erosion  losses  and  the  necessary  repair 
time  are  determined  in  the  EROSIN  routine.  Traffic  delay  losses  are  then 
calculated  in  the  DETOUR  routine.  Following  the  DETOUR  routine,  the 
SCOUR  routine  is  called  to  determine  the  scour  losses  for  the  given  flow 
with  the  backwater  property  damage  determined  in  the  BACKWR  routine.  All 
of  the  above  losses  are  then  summed  and  weighted  and  the  OUTPUT  routine 
is  called  to  print  the  results  for  each  flow.  After  all  flows  have  been 
considered  the  program  exits  the  04  loop.  All  weighted  losses  are  then 
summed  to  give  the  total  risk  for  the  considered  bridge  configuration. 
The  resulting  total  risk  is  printed  in  the  0UTPUT1  routine.  After  all 
bridge  configurations  have  been  considered,  the  a2  loop  is  exited  and 
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a  consolidated  summary  of  results  for  all  bridge  designs  is  printed  in  the 
CONSDA  routine.  The  user  has  the  option  to  plot  the  summary  results  in  the 
PLOT  routine.  Upon  completion  of  the  PLOT  routine  the  program  is  terminated, 


SUBROUTINE  INRISK 

INRISK  is  a  subroutine  called  by  the  main  program  RISK.  This  sub- 
routine reads  and  echoes  all  user  input.  The  necessary  user  input  is  shown 
in  the  flowchart  in  Figure  25.   The  detailed  input  card  formats  are  dis- 
cussed in  subsequent  sections. 


FUNCTION  CONCST 

This  function  evaluates  the  total  cost  of  the  particular  bridge 
superstructure  and  substructure  by  means  of  a  regression  equation  or  a 
user  specified  function.  The  embankment  fill  volume  is  calculated  for 
each  embankment  section.  The  fill  volumes  are  then  added  to  obtain  a  total 
embankment  fill  volume.  This  total  fill  volume  is  then  converted  to 
a  dollar  value  on  the  basis  of  a  per  unit  fill  volume  cost  for  embankment 
fill.  The  length  of  roadway  is  calculated  for  both  embankments  and  over 
the  bridge.  The  roadway  cost  is  determined  by  multiplying  the  roadway 
length  by  a  per  unit  length  cost  factor.  The  bridge  costs,  the  embankment 
fill  costs,  and  the  roadway  costs  are  summed  and  converted  to  an  equivalent 
annual  cost  using  a  specified  project  life  and  interest  rate  by  the 
capital  recovery  factor.  If  the  bridge  has  piles  which  are  driven  to 
bedrock  to  prevent  scour  damage  (TYP=0) ,  this  additional  cost  is  added 
to  the  bridge  construction  cost  computed  above.  The  flowchart  for  this 
function  is  shown  in  Figure  26. 
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(            ENTRY      \ 
\   INRISK/INRISK  j 

v 

© 

I       READ  TITLE  AMD     / 
/HYDRAULIC  INFORMATION/ 

V 

/    DEVICE  NUMBER     / 

/  READ  / 
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f 

^ 

^ 

/     READ  BRIDGE      / 

f 

/      COST  DATA      / 

/  READ  / 
/  EROSION  DATA  / 

f 

j                        DATA/ 

1 

^ 

/        READ          / 
/    EMBANKMENT  DATA    / 

r 

/  PRINT  ALL  / 
/   INPUT  DATA  / 

9 

1 

s 

/       READ         / 
/    BACKWATER  DATA     / 

r 

f     RETURN   ") 

^ 

f 

/       READ         / 
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1 

r 
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^ 

r 
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Figure    25.     Subroutine  INRISK  Flowchart 
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LOOP  TO  l^ 
FOR  EACH  FILL 
SECTION 


ENTRY 
CONCST/CONCST 


CALCULATE 
BRIDGE  COSTS 


CALCULATE 
EMBANKMENT 
FILL  COST 


CALCULATE 
ROADWAY  COST 


SUM  COSTS 


NO 


YES 


ADD  PILING  COSTS 

TO  BRIDGE 

CONSTRUCTION  COST 


RETURN 


CONVERT  TO 
ANNUAL  COSTS 


& 


Figure  26.  Construction  Cost  Subroutine 
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SUBROUTINE  RDHYDR 

This  subroutine  accesses  the  hydraulic  tape  for  each  flow-embankment 
opening  combination.  An  indicator  (ICK)  to  represent  whether  overtopping 
of  the  embankments  has  occurred  is  read  along  with  the  number  of  possible 
overtopping  nodes,  NP.  The  x  and  z  velocity  components  of  flow,  water 
surface  depth,  and  water  surface  elevation  are  then  read  from  the  hydraulic 
tape  for  each  overtopping  node  in  the  particular  configuration.  The  flow- 
chart for  this  subroutine  is  shown  in  Figure  27. 


SUBROUTINE  OVTME 

Subroutine  OVTME  calculates  the  amount  of  time  an  embankment 

section  was  overtopped.  It  then  examines  the  overtopping  time  for  each 

embankment  section  and  the  bridge  itself  to  determine  the  maximum  over- 
topping time.  The  flowchart  is  shown  in  Figure  28  . 


SUBROUTINE  EROSIN 

The  EROSIN  subroutine  calculates  the  amount  of  erosion  by  embank- 
ment section.  It  also  calculates  the  roadway  loss  in  each  section  of  the 
embankment.  It  determines  a  total  percentage  of  fill  volume  lost  by 
erosion  and  uses  this  percentage  to  determine  the  time  to  repair  the 
embankments.  The  losses  due  to  erosion,  namely,  fill  volume  losses  and 
roadway  losses,  are  summed  and  multiplied  by  a  weighting  factor  to  deter- 
mine the  erosion  risk  component.  The  flowchart  for  this  subroutine  is 
shown  in  Figure  29  . 
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ENTRY 
RDHYDR/RDHYDR 


READ  (ICK) 

OVERTOPPING 
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READ  (NP) 
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READ  X,  Z 
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OVERTOPPING  NODE 


RETURN 


Figure  27  .  Subroutine  RDHYDR  Flowchart 
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CALCULATE 

OVERTOPPING 

TIME  FOR  EACH 

SECTION 


DETERMINE 

MAXIMUM 

OVERTOPPING 

TIME 


RETURN 


NO 


SET  OVERTOPPING 
TIME  TO  ZERO 


SET  MAXIMUM 
OVERTOPPING  TIME 
TO  ZERO 


Figure  28.  Subroutine  OVTME  Flowchart 
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LOOP  TO  al  FOR 

EACH  EMBANKMENT 
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ENTRY 
EROS  IN/EROS  IN 


NO 


CALCULATE 

EROSION  LOSS 

OF  SECTION 


SET  EROSION 

LOSSES  AND  DELAY 

TIME  TO  ZERO 
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ROADWAY  LOSS 

OF  SECTION 


DETERMINE  TOTAL 
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COMPONENT 


RETURN 


Figure  29  .  Subroutine  EROSIN  Flowchart 
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FUNCTION  DETOUR 

The  function  DETOUR  calculates  the  cost  of  traffic  delay  caused 
by  embankment  overtopping.  Traffic-related  losses  are  divided  into  four 
categories:  1)  detour  travel  cost,  2)  detour  accident  costs,  3)  unexpected 
obstacle  cost,  and  4)  cost  of  time  delay.  All  four  traffic-related  losses 
are  summed  and  multiplied  by  a  weighting  factor  to  determine  the  DETOUR  risk 
component.  The  flowchart  for  this  function  is  shown  in  Figure  ;30* 


FUNCTION  SCOUR 

This  function  computes  the  pier  scour  damage  caused  by  flooding  of  each 
bridge  and  embankment  configuration.  Three  types  of  bridge  foundations  are 
considered: 

1.  Piles  are  driven  to  bedrock  and  no  scour  damage  is  assumed 
to  occur.  The  data  input  refers  to  this  type  as  Type  "0". 
The  cost  of  piles  should  be  added  to  the  bridge  construction 
cost. 


2.  The  second  type  has  either  the  footing  on  bedrock  or  counter- 
measures  being  taken  to  prevent  pier  scour.  Again  pier  scour 
damages  are  assumed  negligible.  Data  input  refers  to  this 
type  as  Type  "1". 

3.  The  third  type  considers  piles  driven  into  alluvial  soil  where 
the  possibility  of  scour  is  a  certainty.  Data  input  refers  to 
this  type  as  Type  "2".  Maximum  scour,  taken  from  input  data, 

is  compared  with  the  pier  scour  depth  calculated  by  the  function 
SCOUR.  This  calculation  is  derived  from  Laursens  "Scour  at 
Bridge  Crossings", (15). 

If  the  scour  depth  is  greater  than  the  maximum  scour  depth  the  bridge 
is  assumed  lost,  if  not  a  scour  repair  cost  is  determined*  All  scour  costs 
incurred  are  multiplied  by  a  weighting  factor  to  determine  the  scour  risk  com- 
ponent. The  flowchart  of  the  routes  is  shown  in  Figure  31   This  function 
also  estimates  bridge  damage  due  to  submergence  and  includes  it  as  a  component 
of  the  initial  scour  damage. 
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Figure  30.  Function  DETOUR  Flowchart 
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Figure  31.  Function  SCOUR  Flowchart 
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FUNCTION  BACKWR 

This  function  returns  the  backwater  damage  due  to  a  particular 

flood  event  for  a  particular  bridge  configuration.  The  actual  backwater 

damage  calculations  are  performed  manually  by-sthe  analyst  using  the  normal 
runs  of  the  FEM  model  and  consulting  the  corresponding  runs  for  each 

particular  bridge  configuration.  The  results  of  these  manual  calculations 

are  input  to  subroutine  INRISK  which  passes  the  information  to  BACKWR. 

BACKWR  selects  the  appropriate  backwater  damage  and  multiplies  the  damage 

by  a  weighting  factor  and  returns  the  result  to  the  main  RISK  program 

as  the  backwater  risk  component.  The  flowchart  is  shown  in  Figure  32  . 


SUBROUTINE  OUTPUT 

This  subroutine  prints  for  the  given  bridge  configuration  the 
overtopping  time  and  the  time  duration  of  the  resulting  delay  for  the 
given  storm  event.  The  detour,  scour,  erosion,  and  backwater  risk  compo- 
nents and  the  total  risk  component  associated  with  the  storm  are  then 
printed.  The  flowchart  for  this  subroutine  is  shown  in  Figure  33  • 


SUBROUTINE  OUTPUT! 

This  subroutine  prints  for  the  given  bridge  configuration  the 
results  of  the  risk  analysis.  The  total  economic  response,  the  construc- 
tion costs  and  the  total  risk  are  printed.  The  flowchart  is  shown  in 
Figure   34. 


SUBROUTINE  CONSDA 

This  subroutine  prints  a  consolidated  table  of  the  risk  analysis 
results  for  each  bridge  configuration.  The  decision  variables,  the  total 
economic  response,  construction  cost  and  total  risk  are  output  for  each 
bridge  design.  The  flow  and  return  frequency  for  all  storms  used  in  the 


140 


ENTRY 
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Figure   32.   Function  BACKWR  Flowchart 
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ENTRY 
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Figure  33  .  Subroutine  OUTPUT  Flowchart 
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ENTRY 
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RETURN 
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Figure   34.  Subroutine  0UTPT1  Flowchart 
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risk  analysis  are  then  printed.  The  project  life  and  interest  rate  for 
calculating  the  equivalent  annual  cost  of  the  bridge  are  output  next. 
This  subroutine  is  shown  in  the  flowchart  in  Figure  35  . 


SUBROUTINE  PLOT 

The  PLOT  subroutine  provides  the  model  users  with  the  option  to 
obtain  plots  of  total  economic  response  versus  embankment  height  for  a 
fixed  bridge  length,  and  total  economic  response  versus  bridge  length  for 
a  fixed  embankment  height.  The  user  may  bypass  the  plot  option  if  desired, 
The  flowchart  for  this  subroutine  is  shown  in  Figure  36  . 


USER  INPUT  DATA  DESCRIPTION 

The  following  description  indicates  the  exact  FORTRAN  card  format 
to  be  used  by  the  RISK  module.  All  user-supplied  input  is  accessed  on 
logical  unit  5.  In  the  preparation  of  card  input  data  all  integer 
variables  should  be  right  justified  in  the  card  columns  indicated. 
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Figure  35  .  Subroutine  CONSDA  Flowchart 
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Figure   36.  Subroutine  PLOT  Flowchart 
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Card  Sequence: 
Number  of  Cards 

Card 
Columns 

FORTRAN 
Name 

Format 
Type 

Description  of 
Input  Value 

1:1 

5-80 

TITLEC1 ,20) 

A 

Header  Information 

1:2 

5-80 

TITLE(2,20) 

A 

Header  Information 

2:1 

21-25 

ITYPE 

I 

Input  Tape  File 

2:1 

26-30 

IBC 

I 

Bridge  Cost  Function: 
1--User  Supplied  Bridge 

Cost  Function 
0— Model  Supplied  Bridge 

Cost  Function 

3:1 

1-10 

NLEN 

Number  of  Bridge  Lengths 

3:2 

1-10 

FLEN(l) 

Bridge  Length  1 

3:2 

11-20 

FLEN(2) 

Bridge  Length  2 

3:2 

21-30 

FLENC3) 

Bridge  Length  3 

3:2 

31-40 

FLEN(4) 

Bridge  Length  4 

3:2 

41-50 

FLEN(5) 

Bridge  Length  5 

4:1* 

11-20 

B(l.I) 

Coefficient  for  Bridge 
Superstructure  Cost 

4:2 

11-20 

B(2,I) 

E 

Coefficient  2  for  Bridge 
Superstructure  Cost 

4:3 

11-20 

B(3,I) 

E 

Pier,  Footing  and  Abut- 
ment Costs 

4:4 

11-20 

B(4,I) 

E 

Spur  Dike  Cost 

5:1 

11-20 

RWCST 

F 

Dollar  Cost/Foot  for 
Paved  Roadway 

5:1 

21-30 

CY 

F 

Dollar  Cost/Yard3  of  Fill 

5:1 

31-40 

VSLP 

F 

Upstream  Embankment 
Slope  in  % 

5:1 

41-50 

VSLP 

F 

Downstream  Embankment 
Slope  in  % 

5:1 

51-60 

BWID 

F 

Average  Bridge  Width  (ft) 

6:1 

1-10 

NFLFR 

I 

Number  of  Flows 

6:2 

1-10 

FLFR(1,1) 

F 

Flow  1 

6:2 

11-20 

FLFR(1,2) 

F 

Weight  1 

*Repeated  for  each  candidate  bridge  design  for 
if  IBC=0. 


IBC-1 .  Ignore  card  group  3 
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Card  Sequence: 
Number  of  Cards 

Card 
Columns 

FORTRAN 
Name 

Format 
Type 

Description  of 
Input  Values 

6:2 

21-30 

FLFR(2,1) 

F 

Flow  2 

6:2 

31-40 

FLFR(2S2) 

F 

Weight  2 

6:2 

41-50 

FLFR(3,1) 

F 

Flow  3 

6:2 

51-60 

FLFR(3,2) 

F 

Weight  3 

6:2 

61-70 

FLFR(4,1) 

F 

Flow  4 

6:2 

71-80 

FLFR(4,2) 

F 

Weight  4 

6:3 

1-10 

FLFR(5,1) 

F 

Flow  5 

6:3 

11-20 

FLFR(5,2) 

F 

Weight  5 

6:3 

21-30 

FLFR(6,1) 

F 

Flow  6 

6:3 

31-40 

FLFR(6,2) 

F 

Weight  6 

7:1 

1-5 

NMB(4) 

I 

Number  of  Different 
Embankments 

8:1 

11-20 

EMB(4,1) 

F 

Delta  Embankment  Height  1 

8:2 

21-30 

EMB(4,2) 

F 

Delta  Embankment  Height  2 

8:3      / 

31-40 

EMB(4,3) 

F 

Delta  Embankment  Height  3 

9:1 

1-5 

NC(4) 

I 

Bridge  Clearance 

9:1 

6-15 

CEN(4) 

F 

Number  of  Nodes  to  Center 
of  Bridge 

10:1 

1-5 

NNDS(4) 

I 

Number  of  Overtopping 
Nodes 

* 
11:1 

5-10 

NVP(4,40) 

I 

Overtopping  Nodes 

11:1 

11-20 

BTEL(4,40) 

F 

Bottom  Embankment 
Elevations 

11:1 

21-30 

BTL(4,40) 

F 

Distance  Between  Embank- 
ment Nodes 

11:1 

31-40 

WELEV(4,40) 

F 

Top  Embankment  Elevations 

Card  Groups  7  through  11  are  repeated  for  each  bridge  length 
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Card  Sequence:     Card      FORTRAN    Format 
Number  of  Cards    Columns     Name      Type 


Description  of  Input  Values 


12:60, 


max 


21-40 


DAM (60) 


Backwater  Damages 


Note:  Expected  order  is  by  bridge  opening,  embankment  height  and  flow. 


Openi 

ng 

Emban 

kmen 

t  Ht 

Flow 

1 

1 

1 

1 

1 

2 

1 

2 

1 

1 

2 

2 

2 

1 

1 

2 

1 

2 

2 

2 

1 

2 

2 

2 

Thus  for  2  flows,  2  openings,  and  2  embankment  heights,  8  damage  values 
would  be  expected  in  the  above  order.  Damages  for  the  nonconstricted 
case  (no  bridge  in  place)  are  listed  first  by  flow. 


12:61 


1-3 


END 


A 


13:1 

11-20 

PLIFE 

F 

13:1 

21-30 

RATE 

F 

13:1 

31-40 

TC 

V 

14:1 

21-30 

TRF(l) 

F 

14:2 

21-30 

TRF(2) 

F 

14:3 

21-30 

TRF(3) 

F 

14:4 

21-30 

TRF(4) 

F 

14:5 

21-30 

TRF(5) 

F 

Indicates  End  of  Damage 
Information 

Project  Life  (yr) 

Interest  Rate  {%) 

Value  of  Time  ($/hr) 

Average  Traffic  Level 
(Vehicle/day) 

Passenger  Car  Fraction  of 
Total  Traffic  Level  (0-1) 

Commercial  Delivery 
Vehicle  Fraction  of  Total 
Traffic  Level  (0-1) 

Single  Unit  Truck 
Fraction  of  Total  Traffic 
Level  (0-1) 

Gasoline  Truck  Fraction 
of  Total  Traffic  Level 
(0-1) 
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Card  Sequence:     Card      FORTRAN    Format 
Number  of  Cards    Columns     Name      Type 


Description  of  Input  Values 


14:6 


14:7 


15:1 


15:1 


21-30 


21-30 


41-50 


51-60 


TRF(6) 


TRF(7) 


14:8 

21-30 

TRF(8) 

F 

15:1 

1-10 

DRATN 

F 

15:1 

11-20 

PERINN 

F 

15:1 

21-30 

PR0DAN 

F 

15:1 

31-40 

PERINU 

F 

PRODAU 


DRATU 


15:1 

61-70 

DEATH 

F 

15:1 

71-80 

CPERI 

F 

15:2 

1-10 

CPRODM 

F 

16:1 

21-30 

WCL 

F 

16:1 

31-40 

WCH 

F 

17:1 

1-10 

TYP 

F 

17:1 


11-20 


CSCR 


Diesel  Semitrailers 
Fraction  of  Total  Traffic 
Level  (0-1) 

Normal  Travel  Distance 
(Miles) 

Normal  Route  Speed  (mph) 

Normal  Death  Rate 
(Deaths/100  Million 
Vehicle  Miles) 

Number  of  Personal 
Injuries/Death-Normal  Case 

Number  of  Property 
Damages/Death-Normal  Case 

Number  of  Personal 
Injuries/Death-Unexpected 
Obstacle  Case 

Number  of  Property 
Damages/Death-Unexpected 
Obstacle  Case 

Unexpected  Obstacle  Death 
Rate  Multiplier  (Deaths 
Unexpected  Obstacle/Normal 
Rate) 

Cost  of  Death 

Cost  of  a  Personal  Injury 

Cost  of  Property  Damage 

Length  of  Longest  Water- 
course (Miles) 

Elevation  Difference  of 
Longest  Watercourse  (ft) 

0  -  Piles  Rest  on  Solid 

Rock 

1  -  Footing  on  Solid  Rock 

2  -  Solid  Rock  not  Found 

$Cost/ft  of  Pile  for  TYP-0 
$Cost/ft  of  Scour  for  TYP-2 
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Card  Sequence: 
Number  of  Cards 

Card 
Column 

FORTRAN 
Name 

Format 
Type 

Description  of  Input  Value 

17:1 

21-30 

HSCR 

F 

Maximum  Scour  Depth  for 
TYP-2 

17:1 

31-40 

NSCR 

I 

Node  at  Which  Velocity  is 
Used  for  Scour  Calculations 

17:1 

41-50 

WSCR 

F 

Width  of  Rectangular  Pier 
Facing  Upstream 

17:1 

51-60 

NPILE 

I 

Number  of  Piles  for  TYP-0 

17:1 

61-70 

DPILE 

F 

Average  Depth  of  Piles 
for  TYP-0 

17:1 

71-80 

CETA 

F 

Loss  per  Foot  of  Bridge 
Submergence 

18:1 

1-10 

VINCP 

F 

Incipient  Velocity  for 
Erosion  (ft/sec) 

18:1 

11-20 

WT 

F 

Specific  Weight  of  Embank- 
ment (lbs/ft3) 

18:1 

21-30 

FAC 

F 

Cost  Adjustment  Factor 

18:1 

31-40 

TTIP 

F 

Inspection  Time  (hrs) 

19:1 

1-4 

XAXIS 

F 

XAXIS=9999  No.  of  Plots 

Desired 
X=Coordinate  of  Origin 

19:2 

7-10 

YAXIS 

F 

Y=Coordinate  of  Origin 

19:3 

13-14 

ITEL 

I 

Always  1 

19:4 

14-16 

IJOIN 

I 

1-Join  Points  on  Plot 
O-Don't  Join  Points  on 
Plot 

* 
20:1 

1-40 

TOP 

A 

Top  of  Plot  Information 

20:2 

1-40 

BOTTOM 

A 

Bottom  of  Plot  Information 

20:3 

1-50 

SIDE 

A 

Side  of  Plot  Information 

** 
21:1 

1-40 

TOP 

A 

Top  of  Plot  Information 

21:2 

1-40 

BOTTOM 

A 

Bottom  of  Plot  Information 

21:3 

1-50 

SIDE 

A 

Side  of  Plot  Information 

*Include  card  group  20  for  each  bridge  length. 
**Include  card  group  21  for  each  embankment  height. 
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A  listing  of  the  trial  data  set  for  the  example  problem  discussed 
in  Chapter  IV  of  this  volume  is  shown  in  Table  25  as  an  aid  to  the  model 
user  in  formulating  his  own  input  data  set. 


DEFINITION  OF  PROGRAM  VARIABLES 

The  following  variables  are  passed  from  the  FEM  module  to  the 
RISK  module  by  means  of  the  hydraulic  tape. 


VEL 
(4,20) 


ICK         Overtopping  indicator  (1  -  Overtopping, 
2  -  No  overtopping). 

NP  Number  of  overtopping  nodes  in  the  system. 

VEL(1 ,20)     X  =  Velocity  of  flow  at  all  overtopping 

nodes  in  the  system. 

VEL(2,20)     Z  =  Velocity  of  flow  at  all  overtopping 

nodes  in  the  system. 

VEL(3,20)     Depth  of  all  overtopping  nodes  in  the 
system. 

VEL(4,20)     Water  surface  elevation  at  all  nodes 
in  the  system. 


The  following  variables  are  used  in  the  RISK  module  alone: 


NI 

Input  file  number. 

NO 

Output  file  number. 

ITYPE 

Tape  number  of  the  hydraulic  tape. 

IBC 

Bridge  cost  function  option  (1  -  User 
Supplied,  0  -  Model  Supplied) 

TITLE(2,20) 

Heading  information. 
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B(4,4) 


FLFR 
(6,2) 


B(*,l)  Intercept  for  bridge  superstructure 
cost. 

B(*,2)  Slope  for  bridge  superstructure  cost. 

B(*,3)  Pier,  footing  and  abutment  cost. 

B(*,4)  Spur  dike  cost. 

CY  Embankment  fill  cost  ($/cu  yd). 

RWCST  Roadway  cost  ($/ft). 

BWID  Bridge  width  (ft). 

USLP  Upstream  embankment  slope  [%). 

DSLP  Downstream  embankment  slope  (%). 

NC(4)  Number  of  nodes  to  the  bridge. 

CEN(4)  Bridge  clearance  (ft). 

NLEN  Number  of  bridge  lengths. 

FLEN(5)  Bridge  lengths  (ft). 

FLFR(6,1)  Flow  (cfs). 

FLFR(6,2)  Return  frequency  weights. 

NRLFR  Number  of  flows. 

NEMB(4)  Number  of  embankment  heights  for  each 
bridge  length. 

EMB(4,3)  Delta  embankment  heights  for  each  bridge 
length. 

NNDS(4)  Number  of  overtopped  nodes  (overtopped 
nodes  are  also  embankment  section  nodes) 

NUP(4,40)  Overtopping  and  embankment  nodes. 

BTEL(4,40)  Bottom  of  embankment  elevations. 

BTL(4,40)  Distance  between  embankment  and  overtopp- 
ing nodes. 

WELEV(4,40)  Top  of  embankment  elevations. 

NRISK  Number  of  backwater  damages. 

DAM(60)  Backwater  damages. 

PLIFE  Projected  life  {yr) . 

RATE  Interest  rate. 

TC  Value  of  time  ($/hr). 


Repeated  for  each  bridge  design  to  be  considered 
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I 


TRF(ll) 


DRATN 

PERINN 
PRODAN 
DRATU 

PRODAU 

PERINU 

DEATH 

CPERI 

CPRODM 

TRF(1 

TRF(2 

TRF(3 

TRF(4 

TRF(5 

TRF(6 

TRF(7 

TRF(8 

TRF(9 

TRF(IO) 

TRF(ll) 

WCL 
TYP 

CSCR 
CSCR 


Number  of  deaths/100  x  10  vehicle 
mi  1 es . 

Number  of  personal  injuries/death. 

Number  of  property  damages/death. 

Death  rate  multiplier  for  obstacles 
(death  obstacles/deaths  normally). 

Number  of  property  damages,  obstacle 
case/death. 

Number  of  personal  injuries,  obstacle 
case/death. 

Average  cost  of  death  by  accident. 

Average  cost  of  personal  injury. 

Average  cost  of  property  damage. 

Average  traffic  level  (vehicles/day). 

Passenger  car  traffic  level  fraction 
(0-1). 

Commercial  delivery  vehicle  traffic 
level  fraction  (0-1). 

Single  unit  truck  traffic  level  fraction 
(0-1). 

Gasoline  truck  traffic  level  fraction 
(0-1). 

Diesel  semitrailer  traffic  level 
fraction  (0-1). 

Normal  travel  distance  (miles). 

Normal  route  speed  (mi/hr). 

Detour  travel  distance  (miles). 

Detour  travel  speed  (mi/hr). 

Average  occupancy  rate  ($/persons/ 
vehicle) . 

Length  of  longest  watercourse  (miles). 

Type  of  scour  (0  -  Piles  rest  on  solid 
rock;  1  -  Footing  on  solid  rock,  no 
scour;  2  -  Solid  rock  not  found). 

Cost  per  foot  of  pier  for  TYP-0. 

Cost  per  foot  of  scour  for  TYP-2. 
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HSCR         Maximum  allowable  scour  depth  for  TYP-2. 

NSCR         Node  at  which  velocity  is  used  for  scour 
calculations. 

WSCR  Width  of  rectangular  pier  facing  upstream. 

NPILE  Number  of  piles  for  TYP-0. 

DPILE  Depth  of  piles  for  TYP-0. 

CETA  Damage  ($)  per  foot  of  bridge  submergence. 

VINCP  Incipient  velocity  for  erosion  (ft/sec). 

WT  Specific  weight  (lb/cu  ft)  of  embankment 

soil . 

FAC  Cost  adjustment  factor. 

TTIP         Inspection  time  (hr). 
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PROGRAM  LISTING 

The  FORTRAN  IV  program  for  the  RISK  module  follows.  The  program 
was  run  on  the  CDC  6400  unit  at  the  Naval  Ship  Research  and  Development 
Center  computer  facilities,  Carderock,  Maryland. 


10 


15 


20 


25 


30 


35 


PRuGkam  N15*  (INPUT, OUTPUT, TAPfc5«INI»UT,TAi>66«0UT»»UT» 

i   uptn 

CUM«UN/t)L'*TK/HEN(5),NUEN»FUfNl6,2).NFL'H»feNa(<(#J)#NM8(<l}, 
LNNu3(«),*ELEV(U,uO)»NUPU,UOj 
CCKMON/JCOUNT/IDAM 
DIMENSION  TIME(UO) 
DIMENSION  N(},3) 
DIMENSION  3(3.3) 
CALL  INfilSK 
JOAM»0 

DO  250  Il8£»NLtN 
LlMeNMD(U-l) 
DO  2M0  I2«l.Ll" 
CNST««CUNC3T(I1,I2) 
CObl«*U# 

DU  230  I3«  l.NFLFR 
IO*MpID*Mti 
C*tL  KDHYDH(H) 
C   F1*0  Tint  OF  yVEHlCfHING 

CALL  uvfME(U,iz,I3#x»n,TlM6) 

C  KINS)        EKUS10N 

"    CALL    tk05lN(H,I2#U,TlMeiE«»iT«f»J 
C  C*LL     TrlAFFIC    DETUUK 

DTK»0t  1UU«<(1  1,12, 13.X AM,T«P> 
G      FIND    COST'OF    B«1DGE    3COtjR 

8Ch*SCUUSCI1,I2,I3»Cn8TR) 
C   ^IND  COST  OF  BACKWATER  DAMAGE 

bKhS0ACKHH(!i,i2,!3) 

CALL  iuTPUTCll,12,13,XAM,TR»,CUSTR,DTR,8C»6R»»BKrt) 
230  CONTINUE 

CALL  L'UTPT1(11,I2,C08TR,CN8TR.H,8) 
200  CONTINUE 
250  CONTINUE 

CALL  CUN80A(H,8) 

C*Lt  PLQT(K) 

STOP  11 

END 
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